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ABSTRACT 

Composites used in space structures are exposed to both extremes in 
temperature and applied mechanical loads. Cracks in the matrix form, 
changing the laminate thermoelastic properties. The goal of the present 
investigation is to develop a predictive methodology to quantify micro- 
cracking in general composite laminates under both thermal and mechanical 
loading. This objective is successfully met through a combination of 
analytical modeling and experimental investigation. In the analysis, the 
stress and displacement distributions in the vicinity of a crack are 
determined using a shear lag model. These are incorporated into an energy 
based cracking criterion to determine the favorability of crack formation. A 
progressive damage algorithm allows the inclusion of material softening 
effects and temperature-dependent material properties. The analysis is 
implemented by a computer code which gives predicted crack density and 
degraded laminate properties as functions of any thermomechanical load 
history. Extensive experimenation provides verification of the analysis. 
AS4/3501-6 graphite/epoxy laminates are manufactured with three different 
layups to investigate ply thickness and orientation effects. Thermal 
specimens are cooled to progressively lower temperatures down to -184°C. 
After conditioning the specimens to each temperature, cracks are counted on 
their edges using optical microscopy and in their interiors by sanding to 
incremental depths. Tensile coupons are loaded monotonically to 
progressively higher loads until failure. Cracks are counted on the coupon 
edges after each loading. A data fit to all available results provides input 
parameters for the analysis and shows them to be material properties, 
independent of geometry and loading. Correlation between experiment and 
analysis is generally very good under both thermal and mechanical loading, 
showing the methodology to be a powerful, unified tool. Delayed crack 
initiation observed in a few cases is attributed to a lack of pre-existing flaws 
assumed by the analysis. Some interactions between adjacent ply groups are 
attributed to local stress concentrations. These two effects are not captured 
by the analysis due to its global nature. The analysis is conservative in these 
cases and agrees well with data after the observed onset of cracking. 
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CHAPTER 1 


INTRODUCTION 


Advanced composite materials have gained increasing use in aerospace 
structural applications in recent years. Their properties, which include high 
specific stiffness, low coefficient of thermal expansion (CTE), and high 
specific strength, have made composites extremely attractive for 
dimensionally critical space structures such as antennae and waveguide 
supports, solar reflectors, and truss tubes. A lightweight space telescope 
support structure, for instance, may be made very stiff and with a near-zero 
GTE, which would not even be possible with conventional monolithic 
materials. 

1.1 DAMAGE CAUSED BY THE SPACE ENVIRONMENT 

Composites utilized in space structures are exposed to mechanical 
loads and extremes in temperature. The individual truss tubes of one design 
for the NASA Space Station, for example, must be capable of withstanding 
axial loads estimated at ±5.3 kN (±1200 lb) [1]. Moreover, as the spacecraft 
orbits the earth, it moves in and out of the earth's shadow; as a result, a truss 
structure is exposed to temperature changes of up to ±139°C (±250°F) in 
geosynchronous earth orbit [2]. 

When a single composite ply undergoes a change in temperature, it 
expands or contracts in different directions according to its CTE. Rotating 
the ply changes the value of its CTE. If the ply is independent and 
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unconstrained, the expansion or contraction occurs freely, and no ply stresses 
result regardless of ply orientation. Figure 1.1(a) illustrates this 
phenomenon. However, when plies are rotated to different orientations and 
laminated together as in Figure 1.1(b), each ply will not be able to expand or 
contract according to its own CTE due to the constraint of adjacent plies. 
This creates high stresses within the plies. Since the matrix has a lower in 
situ failure strain than the fiber, cracks in the matrix form, called 
microcracks. 

Similarly, when a mechanical load is applied to a composite laminate, 
some of the plies may be loaded in the direction of their fibers, where they are 
strongest. Laminate constraint will cause others to be loaded in their weak, 
matrix- dominated direction, even though those plies may not be bearing 
much of the total applied stress. As a result, microcracks will form in the 
matrices of the latter, as shown in Figure 1.1(c). For both thermal and 
mechanical loading, we define transverse microcracks as cracks which extend 
parallel to the fibers and perpendicular to the x-y plane. Typical microcracks 
are shown in Figure 1.2. 

Microcracking in composites generally does not lead directly to failure. 
The fibers and adjacent plies serve as obstructions to crack propagation, 
preventing a dominant crack from forming as in monolithic materials. It can, 
however, facilitate other damage modes, such as delamination, which could in 
turn lead to premature failure. More importantly, microcracking may cause 
profound property changes, specifically in CTE, stiffness, and Poisson's ratio. 
The structure will subsequently respond differently to future loads. This can 
be detrimental in dimen-sionally critical applications. The design of the 
Space Station truss tubes is driven largely by a CTE requirement of 0 ± 0.9 x 
10'6/°C, and column buckling has been identified as the dominant mode of 
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Figure 1.1 (a) Thermally loaded independent, unconstrained plies; no 

cracks form, (b) Thermally loaded constrained laminate; 
cracks form due to internal stresses, (c) Mec h a ni cally loaded 
laminate; cracks form in transverse ply. 
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Figure 1.2 


Photomicrograph of the edge of an AS4/3501-6 [0^60^60.]. 
laminate under lOOx magnification. Transverse microcrac 
visible in all plies after mechanical loading. 
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failure [1]. Property changes due to microcracking, CTE and stiffness in this 
example, could cause the structure to stray from design criteria, with 
potentially disastrous consequences. 

1.2 PRESENT WORK 

Researchers recognized the importance of microcracking in composite 
structures decades ago. Since then, numerous experimental studies have 
confirmed the resulting property degradation and have identified some key 
factors in microcracking. These include material type, ply thickness, and 
laminate geometry. Even though an extensive data base has emerged from 
this work, it would be unreasonable to conduct tests on every new material 
that is developed and for the infinite number of laminate geometries that are 
possible. Instead, we need a methodology to predict the initiation and extent 
of microcracking analytically. We can then implement existing tools and data 
bases giving property changes as functions of the number of cracks. This 
would also help us to better understand the problem. 

The goal of the present work is to develop an analytical methodology to 
predict the initiation and extent of microcracking in composite laminates 
under both thermal and mechanical loading. The analysis uses a shear lag 
solution of the laminate stress and displacement distributions in the vicinity 
of a transverse microcrack. This solution is subsequently combined with a 
cracking criterion based on fracture mechanics to determine whether further 
cracking is energetically favorable. An incremental damage model facilitates 
the inclusion of temperature dependent material properties and laminate 
softening due to cracking. The analysis is integrated into a computer code 
which can accommodate any user-defined thermomechanical load history. 
Extensive testing, using a variety of laminate geometries under thermal and 
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mechanical loading, was performed to verify the analysis. The data collection 
included crack counts in the specimen interiors as well as at their edges. 

A data fit provided two key parameters for the analysis, both assumed 
to be material properties. Some interaction between cracks in adjacent plies, 
not included in the analytical model, was observed. In general, though, 
correlation between analysis and experiment was excellent under both 
thermal and mechanical loading and for all laminate geometries. 

1.3 OVERVIEW 

In Chapter 2 the previous work relevant to the problem is reviewed. 
This includes work which first recognized the problem, studies leading to the 
present work, and recent work on the problem. A problem statement 
summarizing the present work is presented in Chapter 3. The analytical 
methodology, including the shear lag model and the formulation of the energy 
expressions, is developed in Chapter 4. Chapter 5 describes the thermal and 
mechanical tests which were used to verify the analysis. The experimental 
results and correlation with predicted values are discussed in Chapter 6. 
This chapter also includes results of the data fits, parametric studies, and 
discussion of strengths and weaknesses of the model. Finally, conclusions 
and recommendations for future work are presented in Chapter 7. 
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CHAPTER 2 


PREVIOUS WORK 


Microcracking has long been recognized as an important damage 
mechanism in composite laminates. Early studies identified the problem and 
described it qualitatively, but most of these concentrated on the effects of 
microcracking rather than on its formation. Later studies developed 
predictive methodologies for both monotonic and cyclic mechanical loading. 
Recently work has been reported on the quantitative prediction of 
microcracking under thermal loading. The vast majority of studies to date 
have concentrated on the specific problem of cross-ply laminates under 
mechanical loading, a small subset of the current problem. No completed 
analytical methodology has demonstrated the ability to quantitatively predict 
microcracking under both thermal and mechanical loading in general angle 
laminates. 

2.1 EARLY MICROCRACKING STUDIES 

One of the first works to examine the implications of microcracking on 
dimensional stability in the space environment was that of Camahort, 
Re nnhack, and Coons [3]. As with most of the early studies, their work was 
essentially phenomenological in nature. Near-zero- CTE composite specimens 
of five different materials were thermally cycled and then measured for CTE. 
Other specimens were exposed to progressive mechanical loading to evaluate 
residual strain as a function of applied stress. Their work confirmed that 
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near-zero-CTE laminates are indeed sensitive to thermal environments. 
They also showed that microcrack formation is highly dependent upon the 
material system. 

Garrett and Bailey [4] studied transverse microcracking in cross-ply 
lamina tes under me chanical loading. They used a one-dimensional shear lag 
model to calculate the stress transferred from the cracked transverse plies to 
the uncracked longitudinal plies. Expressions for the degraded laminate 
properties as functions of the applied load and uncracked properties were 
derived. The theoretical effective moduli just before failure correlated well 
with experimental results. This work played a key role in characterizing the 
microcracking problem and emphasized the need for predictive 
methodologies. 


2.2 PREDICTIVE METHODOLOGIES 

A number of analytical approaches exist to quantitatively predict 
microcracking in composite laminates. All of the various methods first 
require knowledge of the stress distribution. This is typically determined 
analytically from one of the following: shear lag analysis, the self-consistent 
method, variational approaches, continuum mechanics, or classical laminated 
plate theory. The stress state can also be found numerically through finite 
element analyses. 

Next, after the stress distribution is determined, a cracking criterion is 
introduced to predict whether new cracks will form. Failure criteria 
generally follow either the strength of materials approach or the fracture 
mechanics approach. Using a strength of materials criterion, the derived 
stress state is compared with material strengths, such as the transverse 
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failure strength of a unidirectional ply. A fracture mechanics cracking 
criterion, on the other hand, is based upon energy methods. From the stress 
distribution model, the amount of energy dissipated when a hypothetical new 
crack forms may be determined. If this calculated energy is greater than the 
critical strain energy release rate, a material property, the criterion for new 
crack appearance is satisfied. 

Variations of the two main approaches are also used. Statistical 
methods may be incorporated into both. In the fracture mechanics approach, 
the location of a new crack may be assigned a probabilistic distribution or 
determined from a Monte Carlo simulation. With the strength based 
approach, a probabilistic distribution of the failure strength may be assumed. 

Studies have also considered the problem on different scales. For 
example, some studies examine microcrack formation on the microstructural 
level, considering the fiber, matrix, and fiber-matrix properties individually. 
However, transverse microcracking analyses generally focus on the ply level, 
at which the fiber, matrix, and interface are considered homogenized. 

2.2.1 Mechanical Loading Analyses 

In the simplest example of a strength based analytical model, the 
stress state is calculated with CLPT, and the criterion for crack initiation is 
the transverse strength of a single lamina [5-8]. This approach relies on the 
assumption that in situ ply strength is a material property. The validity of 
this method was called into question by Flaggs and Kural [9], who showed 
that in situ lamina transverse strength is instead a laminate property, 
dependent upon laminate geometry, stacking sequence, and ply thickness. 

The strength criterion has also been used by Parvizi and Bailey, along 
with a shear lag solution of the stress state [10]. Lee and Daniel combined 
this fail ure criterion with a "modified” shear lag model, which assumed a 
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more complex through-thickness displacement profile [11, 12]. Allen et al. 
determined the stress fields in a cracked laminate using the internal state 
variable concept, a continuum damage mechanics approach [5, 13]. In their 
method the damage, transverse microcracking in this case, was represented 
by a second order tensor. Calculated stresses were compared with transverse 
ply strength. While qualitatively useful, these approaches are still limited by 
the laminate dependence of the in situ transverse ply strength. 

Peters et al. [14-17] used a shear lag stress solution combined with a 
two-parameter Weibull strength distribution as the cracking criterion. Jen 
and Sun [18] incorporated a similar strength distribution with a finite 
element stress solution. The strength and shape parameters of the Weibull 
distribution are determined experimentally from a single specimen. These 
analyses have limited application in general cases, as the Weibull parameters 
are related to the in situ transverse ply strength and are thus laminate 
dependent. The works are useful parametric studies, however, since the 
effects of laminate geometry, stacking sequence, and ply thickness on the 
Weibull parameters help to better characterize the transverse microcracking 
problem. 

Many researchers have used a fracture mechanics failure criterion 
with a shear lag stress solution to predict transverse microcracking in cross- 
ply laminates [19-26]. Some, such as that of Caslini et al., predict only the 
onset of microcracking [19]. Laws and Dvorak, in contrast, proposed a 
progressive damage model which determines both the onset and 
accumulation of cracks. It also incorporates a probabilistic distribution for 
the site of the next crack [24]. 

Hashin pioneered the variational approach to determine the 
thermoelastic properties of a cracked laminate [27]. In this method the 
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complementary energy is minimized to derive the stress and strain fields. 
Though his analysis cannot predict the onset or accumulation of 
microcracking, Naim [28, 29] and Varna and Berglund [30] later combined 
his stress analysis technique with a fracture mechanics failure criterion to 
predict microcracking in cross-ply laminates. 

Other energy based approaches include studies by Wang and 
Crossman [31] and Binienda et al. [32], who calculated stress distributions 
with finite element models. Nuismer and Tan developed a two-dimensional 
elasticity model for their analysis [33]. Laws and Dvorak [34, 35] proposed a 
self-consistent model in which the cracking layer is replaced by an effective 
homogeneous medium containing many aligned slit cracks. Wang et al. [36] 
used the concept of effective flaw distributions in conjunction with a Monte 
Carlo simulation. The size and spacing of the microflaws are determined by 
probability distributions. A new crack is assumed to form at the most 
significant microflaw. As a result, this method is often called the "weakest 
link" approach. 

In general, the applicability of predictive methodologies for mechanical 
loading is confined to monotonic loading of cross-ply laminates. Only a few 
analyses have been derived for general angle-ply laminates. Gudmundson et 
al. used a variational approach to determine the stress fields and 
thermoelastic properties of a cracked laminate [37]. Flaggs used a shear lag 
stress analysis with a mixed mode fracture mechanics failure criterion to 
predict microcracking initiation [20]. These analyses, however, have not 
demonstrated the ability to predict microcrack accumulation as a function of 
a progressively applied load. 
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2X2 Thermal Loading Analyses 

Very few predictive methodologies exist for thermal loading. A number 
of the mecha n ical loading analyses incorporate a residual thermal stress due 
to manufacture [5, 11, 13, 15, 16, 23, 24, 26-29, 31, 32, 34-38]. None of them 
have demonstrated predictive capabilities for progressive thermal loading. 
The residual stress is generally treated as a constant. 

Thermal loading analyses have generally focused on issues other than 
prediction of cracking, such as microcracking effects. Bowles [39] and Adams 
and Herakovich [40] used finite element stress analyses to determine the 
effects of microcracking on longitudinal CTE. Bowles also proposed a 
technique to use finite element analysis in conjunction with CLPT to 
determine all of the thermoelastic properties, though no experimental 
verification was presented [39]. Tompkins et al. [41-44] and Camahort et al. 
[3] conducted experimental studies in which the effects of thermal cycling on 
laminate properties and microcrack density were measured. Manders and 
Maas [45] tested thin plies and Bowles and Shen [46] thin fabric to determine 
whether microcracking would be suppressed. Experimental testing by Knouff 
showed that the accumulation of microcracks under thermal cycling is 
dependent upon fiber type [47]. 

Some studies have tried to quantitatively predict microcracking under 
thermal loads. Bowles et al. [46, 48] and Adams et aL [49] used CLPT stress 
solutions with in eitu transverse strengths to predict cracking initiation. A 
similar approach was taken by Peters and Andersen [38], though with 
Weibull distributions of the failure strengths. Chamis et al. developed a 
finite element model which calculates nodal stress resultants in complex 
structures under thermal loading. A combined stress failure criterion is used 
to simulate damage progression [50]. Herakovich and Hyer predicted 
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initiation temperature with CLPT and microcrack accumulation with finite 
elements in cross-ply laminates subjected to thermal loading; both stress 
analyses use an in situ ply strength failure criterion [51]. These approaches, 
besides having limited application to microcrack accumulation under 
progressive loading, suffer from the problems associated with strength based 
failure criteria. 

A few analyses have focused on damage at the microstructural level 
[48, 52-55]. Their objective is to determine the effects of constituent 
properties, interface strengths, and other micromechanical parameters on the 
development of microflaws. Using finite elements, Bowles and Griffin [53, 
54] modelled the thermally induced stresses associated with continuous fibers 
embedded in a matrix. The location and initiation temperature of fiber 
debonding is predicted by comparing the radial stress at the fiber-matrix 
interface with the radial tensile strength of the interface. 

2.3 RECENT WORK 

McManus et al . [56] developed a predictive methodology for cross-ply 
laminates under thermal loading. They used a shear lag stress solution with 
an energy based failure criterion. Crack density in one 90° ply group and 
knocked-down laminate properties were derived as functions of progressively 
decreasing temperature. An iterative degradation model included the effects 
of thermal cycling. A computer code was written to implement the analysis, 
giv ing crack density and reduced laminate properties as functions of 
progressive temperature change or number of thermal cycles. 

They also collected experimental data to correlate with the analytical 
predictions. Cross-ply specimens were progressively cooled and thermally 
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cycled in an environmental chamber. Microcracks were counted on the edges 
of the specimens using optical microscopy. 

The analysis correlated well with both the progressive cooling and 
thermal cycling data. One significant finding was that, in thermal cycling, 
the thinnest ply groups cracked much earlier and more extensively than 
predicted. This result is contrary to recent trends to suppress microcracking 
by using thinner plies. However, cracking in thin plies generally has little 
effect on the effective laminate properties. 

This work was extended by Park [57] to predict transverse 
microcracking in all plies of a general angle-ply laminate. The incremental 
damage model was also modified to include the effects of temperature- 
dependent material properties. The existing computer code was changed 
extensively to incorporate the new approach. Crack densities were measured 
experimentally in a variety of laminates exposed to progressive thermal 
loading. Cracks were counted at the specimen edges using optical 
microscopy, and in the interiors using x-ray inspection and by incrementally 
sanding down the edges. 

Edge crack counts in thin plies did not correlate well with the 
analytical predictions. The agreement between theoretical crack densities 
and interior crack counts, though, was significantly better. The poor 
correlation of the edge crack counts was attributed to two effects. First, 
observations from the interior data showed that cracks in the thin ply groups 
were discontinuous, but the model assumes that transverse cracks form edge- 
to-edge. Second, a free- edge stress analysis showed that the stresses are 
lower in this region. Thus fewer microcracks will form at the free edge than 
in the interior. As a result, edge crack counts are probably not a good 
indication of the damage state within the laminate. 
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An important finding of the work was that correlation between 
analysis and the average of the interior crack measurements was very good, 
even though scatter was high. Thus the analysis appears to be a useful tool 
for predicting mean crack densities under thermal loading, even for thin 
plies. 

2.4 SHEAR LAG AND FRACTURE MECHANICS PARAMETERS 

Two parameters are necessary when a shear lag stress solution is used 
in conjunction with an energy based cracking criterion to predict 
microcracking. The first is the critical strain energy release rate, G^, a 
measure of fracture toughness. This parameter is physically meaningful but 
currently impossible to measure directly for microcracking. Many studies 
instead use the interlaminar fracture toughness, which can be measured 
using a double cantilever beam test. However, this property characterizes 
delamination formation and growth, which is an entirely separate problem 
from microcracking, or intralaminar fracture. Moreover, both the accuracy 
and precision of current measuring techniques are questionable. Figure 2.1 
shows the broad range of G /e , values that were measured during round-robin 
testing in one study [58]. 

The second is a non-dimensional shear lag parameter, referred to as £ 
in the present study. Note that a standardized notation does not exist for 
this parameter. Furthermore, comparison of values from different studies 
must be done carefully. For example, some analyses use a shear lag 
parameter with dimensions of length. This semi-empirical aspect of the 
shear lag model is both a material property and a geometric factor. Like G Ic , 
this parameter is impossible to measure directly. Some have assumed the 
shear lag parameter to be laminate dependent. Laws and Dvorak [24], for 
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Figure 2.1 After [68]. Results of round robin G h testing in three different 
laboratories using four different test techniques. 
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example, back-calculated £ from first-ply-failure and interlaminar fracture 
toughness data. Runkle [59] used a similar approach and attributed the 
laminate dependence of £ to local delamination. Lee and Daniel [12] and 
Ogin et al. [60] derived closed-form solutions for £ based upon kinematic 
assumptions. Others have determined £ empirically, assuming its value to be 
fixed for a given fiber type and independent of laminate geometry [15, 16, 56, 
57, 61, 62]. 
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CHAPTER 3 


APPROACH 

3.1 PROBLEM STATEMENT 

The objective of the present investigation is to develop an analytical 
methodology to predict microcracking in composite laminates under both 
thermal and mechanical loading. The motivation for the present work is 
provided by the lack of understanding in the area of microcrack formation as 
identified in Chapter 2. Analytical modeling in conjunction with 
experimental investigation is used to achieve this research objective. 

The problems addressed by the analysis are threefold: (i) to provide a 
Haaig n capability such that given material properties, laminate geometry, and 
thermomechanical load history, the analysis will predict crack densities and 
degraded laminate properties; (ii) to aid in identifying key parameters in the 
microcracking problem as well as the sensitivity to these parameters; and 
(iii) to provide insight to the physical mechanisms involved in crack 
formation. 

The purposes of the experimental investigation are (i) to provide input 
parameters for the analysis; (ii) to provide verification of the analysis; and 
(iii) to allow a greater qualitative understanding of the problem, including the 
effects of width, layer thickness, and laminate geometry. 
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3.2 ANALYTICAL APPROACH 

The analysis is a fracture mechanics approach which uses a shear lag 
derivation of the stress and displacement fields. All plies in a general 
composite laminate are analyzed under progressive thermal and mechanical 
loading. Each ply is analyzed individually with the properties of the other 
layers smeared. The damage is modeled incrementally, allowing the 
inclusion of existing cracks in other layers and temperature dependent 
material properties. 

A transverse microcrack is assumed to form parallel to the fibers in a 
respective layer and to extend completely through both the layer thickness 
and the laminate width. The analysis assumes that a crack appears 
instantaneously when the conditions for crack formation are met. A new 
crack forms when the energy released due to crack formation is greater than 
the energy required to form a new crack surface, the critical strain energy 
release rate G Ie . The crack formation process is illustrated in Figure 3.1. 

The analytical methodology is implemented by a computer code. A 
parametric study is used to determine the sensitivity of the analysis to the 
various input parameters. Full derivation of the analysis is presented in 
Chapter 4. 

3.3 EXPERIMENTAL APPROACH 

All laminates are fabricated at MIT Technology Laboratory for 
Advanced Composites (TELAC). AS4/3501-6 graphite/epoxy is chosen as the 
material system to build on the data base established by the laboratory. 
Three different layups are manufactured to identify the effects of laminate 
geometry and ply thickness. 
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Figure 3.1 (a) Laminate containing transverse microcracks. Cracks run 

completely from edge to edge and through thickness, (b) 
Cracks assumed to appear instantaneously when conditions for 
new crack formation are met. (From Park, [57]). 
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Two different sizes of thermal specimens are chosen to investigate 
width effects. The mechanical specimens are the same width as the wider 
thermal specimens. Prior to testing, both long edges of every specimen are 
polished to minimize crack initiation sites and to facilitate microscopic 
inspection. Thermal specimens are progressively cooled from room 
temperature to -184°C in a thermal environment chamber. After cooling to 
each of five target temperatures, cracks are counted using optical microscopy. 
Some thermal specimens are used for edge crack data. Others are sacrificed 
to collect interior data by sanding down the edges. The mechanical 
specimens are loaded under monotonic tension. Edge crack data is collected 
at progressively higher loads until failure. 

A fit of the data provides shear lag and fracture toughness parameters 
for the analysis. Analytical predictions are subsequently correlated with test 
results to validate the model and provide insight into the mechanics of the 
problem. 
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CHAPTER 4 


ANALYTICAL METHODS 


In this chapter the derivation of the analysis and its implementation 
are presented. The analysis uses a shear lag solution of the stress and 
displacement fields and an energy based cracking criterion. All plies in a 
general composite laminate are analyzed under progressive thermal and 
mechanical loading. Also described is a data fit procedure used to obtain key 
parameters for the analysis from experimental data. 

4.1 DERIVATION 

4.1.1 Shear Lag Stress Solution 

Figure 4.1 shows a laminate aligned with a global coordinate system 
xyz. The laminate is made up of unidirectional plies. Stacked plies with the 
same ply angle are assumed to act as a single thick ply, referred to as a ply 
group or layer. Cracks are assumed to span the ply group thickness and 
propagate parallel to the fibers through the width of the composite laminate. 
Figure 4.1 shows cracks of this type in a laminate with three ply groups. A 
local coordinate system x'y'z ' is defined for each crack. The y ' axis is aligned 
with the crack, parallel to the fiber direction of the ply group, the x axis is 
aligned with the transverse direction of the ply group, and the origin is at the 
center of the crack. 

To predict cracking in any one ply group, the laminate is modeled as 
being made up of two components: the cracking ply group and the rest of the 
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Figure 4.1 Global and local coordinate systems and assumed microcrack 
geometry. 


42 






laminate, which is smeared. A one-dimensional shear lag model with unit 
depth, shown in Figure 4.2, is used to determine the stress and displacement 
distributions in the vicinity of a crack in the local x'y'z' coordinate system. 
The model assumes uniform through-thickness displacement and normal 
stresses in every layer. Shear stress exists only in a shear transfer region of 
thickness a, between layers and is uniform through-thickness within that 
region. 

In the shear lag solution that follows, E,a, u, and <x are the stiffness, 
CTE, displacement, and normal stress, respectively, in the x' direction. The 
subscripts c, r, and o indicate the cracking layer, the rest of the laminate, and 
the entire laminate, respectively. The shear stress between uncracked and 
cracked layers in the x'y' plane is q. The thicknesses of the cracking layer 
and the rest of the laminate are a e and a„ respectively, and a« is the total 
laminate thickness. 

The laminate is subjected to a thermal load, AT, which is the difference 
between the current temperature and some stress-free temperature, and an 
applied stress, <?„ where 



and N , is the laminate load in the x‘ direction. N , is determined by 

X X 

transforming N, the laminate load vector in the global coordinate system 
used in CLPT [63], to the local coordinate system of the cracking layer. From 
equilibrium of the laminate, 

<7 a a 0 = a r a r + a c a e (4.2) 

From equilibrium of the cracked layer. 
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Figure 4.2 Shear lag model used 
distributions. 


o c + (da/dxOdx' °o 

__»<7 


o r + (do/dxOdx' 


to derive stress and displacement 




(4.3) 


<7 = 


2 dx' 


and of the rest of the laminate, 


a, da, 
2 dx' 


(4.4) 


The stress-strain relations for the cracked layer and the rest of the laminate, 
respectively, are 


and 


Let 



du c 

dx' 


-a e A T 



du, 

dx' 


-a r AT 


(4.5) 


(4.6) 


g = K(u e -u r ) (4.7) 

where if is a stiffness constant relating the displacements of the two layers to 
the shear stress between them. It can be determined from 



(4.8) 


where G*^ is the effective shear modulus of the shear transfer region. 
Placing Eq. (4.7) into Eqs. (4.3) and (4.4) and taking the derivative of both 
gives 


and 


\dx' dx ' ) 2(d x'f (4.9) 


Idx' dx' J 2 (dx') a (4.10) 
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Subtracting Eq. (4.6) from Eq. (4.5), multiplying by K, and substituting into 
Eq. (4.9), 


(dx') a e E e E r K rJ J 


(dx 

Solving Eq. (4.2) for o r , substituting into Eq. (4.11), and rearranging. 


d 2<y e o K 

'a r E' + aJl.) 

2K 


a r a e E& ) 

<T,= ‘«r 


Let the shear lag parameter £ be expressed as 


r- I EEsBsEs 

* V 2a r E r E e 


Also let 




It will be convenient to express this parameter as 
where < 7 , , the far-field stress in the cracking layer, is 




Es. 

L*. 




From the rule of mixtures. 


a 0 E t = a^,. + a e E e 

Substituting Eqs. (4.13), (4.14), and (4.17) into Eq. (4.12) gives 

s-x 

(dje ') a <■’ * 

This has a solution of the form: 


(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 
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a. = Asinh — 1 — I + Bcoshl — 1 — + 


a e J 4 C 2 


(4.19) 


Applying boundary conditions cr e = 0 at x'= ±h gives 


A = 0 


4£ a cosh(2£/i/a c ) 


(4.20) 


Placing Eqs. (4.20) into Eq. (4.19), substituting Eq. (4.15) into the result, and 


rearranging, 


L.coshg^l 


(4.21) 


* L cosh(2W« c )f- 

It can be seen by examination of Eq. (4.21) that the shear lag parameter, 
scales the distance from the crack face over which the stress rises towards its 
far-field value. Solving Eq. (4.2) for a r , substituting into Eq. (4.21), and 


rearranging gives 


- - a a a 0 f^Jl C08h( 2 
f a r ~ a r cosh(2CA/o e )J ** 


(4.22) 


Substituting Eq. (4.16) into the equation formed by placing Eq. (4.21) into Eq. 
(4.5), solving for du e , and integrating from 0 to x\ 


u 

e ° 0 E 0 


[o^ J[ 2£x cosh(2Ch/a e ) 


.„ A7V i 5A[i a. sinh(2£r/q c ) 
A1 * a 0 E 0 1 1 2^' cosh(2£A/a e ) 


(4.23) 


Substituting Eq. (4.22) into Eq. (4.6), solving for du r and integrating from 0 


to x\ 


It , P f. a. «jnh(2{»yq.) V| 

r q,£AL * * ,£ V"2C«'co8h(2Ch/q.)JJ 

A7V L f f / / t a, sinh(2Qc7aJ N j1 

a,E. [ • E ’ ' a ‘ E '( a < a ^ 1_ 2ijc' coah(2CV°.) JJ 
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The constants of integration in Eqs. (4.23) and (4.24) are set to zero to satisfy 
the displacement boundary conditions u r ( 0) = « c (0) = 0. 

4.1.2 Energy Expressions 

Two criteria must be met for a new crack to form. First, the crack 
must be energetically favorable. The strain energy release rate associated 
with crack growth, G, must reach a critical value G ICt viz 

G ^ G fc (4.25) 

where 

G = ±±(W-U) (4.26) 

a e 61 

W is the work done by an external load, U is the strain energy stored by the 
body, / is the crack length, and a e is the thickness of the cracking body [64]. 
G Ie is the Mode I critical strain energy release rate, a material property, also 
referred to as the fracture toughness. Only the Mode I fracture toughness is 
considered. The x '-direction stresses calculated in the one-dimensional shear 
lag model are Mode I opening stresses. A similar model could be used to 
calculate the Mode II x'y' shear stresses. However, Mode II fracture is 
considered a secondary effect in the present investigation and is thus not 
included in the analysis. The Mode II fracture toughness, G Uo is generally 
significantly greater than G*, and the x'y' shear stresses are generally lower 
in magnitude than the Mode I opening stresses, especially under thermal 
loading. 

The second criterion for crack formation is that it must be 
mechanistically possible. The strain energy release rate in Eq. (4.26) is a 
function of crack length. Due to the constraint of adjacent plies, G is 
independent of crack length after a critical crack size l« is reached [57], as 
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shown in Figure 4.3. Thus the strain energy release rate may be calculated 
from the total change in energy between uncracked and cracked states in a 
Griffith energy balance [64]: 

G = AW ~* U (4.27) 

<*c 

where G is the strain energy release rate and has units of energy per unit 
area. AW and AU are the changes in external work and internal strain 
energy, respectively, from the state of Figure 3.1a to that of Figure 3.1b. 
Equation (4.27) is derived for the case of unit depth, which is implied in the 
derivation of all energy expressions herein. If the cracking layer is assumed 
to contain initial starter cracks of length l e [64], then Eq. (4.27) may be 
substituted into Eq. (4.25) to get 


AW -ALT 


*G te 


(4.28) 


which is the condition for the appearance of a new crack. 

The change in internal strain energy has contributions from both the 
normal and shear stresses. The strain energy from normal stresses, is 

= (4.29) 

9 

The cracks are initially separated by a distance 2 h. This is shown in the 
volume element in Figure 4.4. The element is aligned with the local x'y'z' 
coor dinat e system of the cracking layer and has emit depth. It will be used in 
the derivation of all energy expressions presented here. A hypothetical new 
crack is assumed to form midway between the existing cracks, such that the 
crack spacing becomes h. The change in U a when a new crack appears is 
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Figure 4.3 Graph of the strain energy as a function of crack length [57]. 
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□ Cracking Layer 
■ Rest of Laminate 


Figure 4.4 Volume element used in energy calculations, showing existing 
crack spacing 2 h. Hypothetical new crack forms midway 
between existing cracks, reducing crack spacing to h. 
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Atf„=[2lU-£/„|J (4 30) 

2 U a \ h is the strain energy due to normal stresses of the volume element with 
cracks separated by a distance h. is the strain energy with cracks 

separated by 2 h and can be found from 


U'L * 2 


|o, f '+±a, f* s±ix' 

2 J. 2 «J„ E, 


ir 


(4.31) 


Note that due to symmetry, we can integrate from 0 to h and multiply the 
result by two, rather than integrating from -h to h. Substituting Eq. (4.22) 
into I, using the trigonometric relation cosh 2 x = ^(cosh2x + l), and 

a 

integrating, 


-a P g 2 
2 'J 0 E r 2a r E, * 
*» 


A-%tanh 


(?) 


-^-sinh( 

8C l 

A Ch' 

s«« * 

2 

cosh 2 

f 2£h 

< a c > 

| 




(4.32) 


Substituting Eq. (4.21) into II, using the trigonometric relation 
cosh a x = ^ (cosh2x + 1), and integrating, 

m 



Placing Eqs. (4.32) and (4.33) into Eq. (4.31), 
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VI .2^ 
*“ a,E r 


+ <V»A 


drEcE, 




/ 

^£-8inh 

8C l 

4CA" 

+ f 

cosh 2 

2£A 

a c . 



<rt 


(4.34) 
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The strain energy due to normal stresses in the volume element after a 
new crack forms and the cracks separated by a distance h is 


[ i i** 3 i f*^ /r* a 

W. t*'+W. t*' 


(4.35) 


Due to symmetry, we integrate from 0 to /i/2 and multiply the result by two. 
Also, in Eq. (4.35), a c and a* are the normal stresses in the cracking layer 
and the rest of the laminate, respectively, derived with the boundary 
conditions a* = 0 at x’ = ±%. Then, by analogy with U^, 

v\ - < * a * h 

" u 2«B 


a cl E 




2o e o. 
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) 





(4.36) 


Substituting Eqs. (4.36) and (4.34) into Eq. (4.30) and using the 
relation smh2r = 2sinhrcoshx, 
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ol (4.37) 




Substituting Eq. (4.16) into Eq. (4.37), expanding the loading terms, and 
factoring, 

'*a)Y 

2 Ca r E r E 0 [ UeJ 2 l <*e J <*c{ l<* e J V <*e )) 

,o r o e 2 S r fi e (q e -g r ) 2 Ar 2 ~ 


a 2 e E c {a e -a r )g a 6 T 


2 CE 0 


2 V a c J l^aj a e ^ l,a e J ^ a e J 

2 tanh|J2 j - tanh 


-sech 2 


(4.38) 


The strain energy contribution from shear stresses, U q> is given by 


(4 - 39) 

v 

Recall that shear stresses are assumed to exist only in the shear transfer 
region between layers. Thus Eq. (4.39) is integrated over the volume of this 
region alone, which has a total thickness of 2a, since shear is transferred at 
both the top and bottom of the cracking layer. The change in U q when a new 
crack forms is 


(4 - 40 > 

2U q \ h is the strain energy due to shear stresses of the volume element with 
cracks separated by a distance h, 17, the strain energy with cracks 
separated by 2 h, can be found by solving Eq. (4.8) for G* and substituting 
the result into Eq. (4.39), giving 
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(4.41) 


U„ 



Due to symmetry, we can once again integrate from 0 to h and multiply the 
result by two, rather than integrating from -h to h. Placing Eqs. (4.23) and 
(4.24) into Eq. (4.7) and squaring, 


Q 


a _ KW a E! 

“ a 2 E 2 E 2 



(4.42) 


Substituting Eq. (4.42) into Eq. (4.41), integrating, and substituting Eq. 
(4.15) into the result. 




Jv 

«C ) 
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AalElEX 

cosh 3 

(Wi 

l a c J 
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(4.43) 


The strain energy due to shear stresses in the volume element when a new 
crack forms and with the cracks now separated by a distance h is 


U a 




(4.44) 


Due to symmetry, we integrate from 0 to hi 2 and multiply the result by two. 
Also, in Eq. (4.35) q*is the shear stress derived with the boundary conditions 
a c ■ 0 at x'=±%. Then, by analogy with , 



KalalEl 
4 a?EX? 




(4.45) 
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Solving Eq. (4.13) for K, placing the result into an equation found by 
substituting Eqs. (4.45) and (4.43) into Eq. (4.40), substituting Eq. (4.16), and 
using the relation sinh2x = 2sinhxcoshx, 


W, = jj^- (a, - a r )6T J 


(4.46) 


The total change in strain energy when a crack appears, AU, is simply 

AU = AU a + AU, (4.47) 

Substituting Eqs. (4.37) and (4.46) into Eq. (4.47) yields 


2tanhf — 

-1-tanh 


U 

c) 

°c J. 


Note that the coupling between thermal and mechanical loading in Eqs. 
(4.37) and (4.46) drops out when the two are added. As a result, the total 
change in internal energy in Eq. (4.48) has no thermomechanical 
interactions, only pure thermal and pure mechanical loading contributions. 

The change in external work when a crack forms, AW, is given by 

AW = [2W| a -WU (4.49) 

where and 2W| Jk are the work done by the applied loading before and 
after the new crack forms, respectively. The former is found from 


n*= 2 w,w 

Setting x'=h in Eq. (4.24) and substituting into Eq. (4.50) gives 


(4.50) 




~ a^.E.aj - a!a r E'E, [a, -a , ) ATcr. 1 2jh tnn( f 2{aV 

1 a c l °c ). 


(4.51) 


+ 2a 0 a r ATaJi 
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The work done by the applied loading after the new crack forms is given by 


n = 2a.a.«;(f) (4.52) 

where is found from Eq. (4.24), derived using the boundary conditions 

< 7 * = 0 at x' = ±%. By analogy with 


n- 


_ ~ a ' a r E ' E '( a ' ~ SL i^-tanhf— ) + a 0 a r ATcr a h 

[ frErE. i°« UJJ 


(4.53) 


Substituting Eqs. (4.53) and (4.51) into Eq. (4.49), 


AW _ 2a e 2 a 0 2 E c c^ - 2a?a, 

2 Co, 


^Vtanhf^ 

'a r a 0 E r E 0 [ a e J a e 


(4.54) 


The change in external work in Eq. (4.54) has thermomechanical coupling 
and pure mechanical loading terms. The amount of work done is the 
displacement u r multiplied by the magnitude of the applied force. If the 
laminate is simultaneously exposed to a change in temperature, u r will have a 
different value than with mechanical loading alone. The work done will 
chang e as well, hence the thermomechanical interaction term. Thermal 
loading alone does no external work and thus has no contributions in Eq. 
(4.54). 

Finally, the strain energy release rate is found by placing Eqs. (4.54) 
and (4.48) into a Griffith energy balance, Eq. (4.27), and substituting Eq. 


(4.16) into the result, yielding 


2tanh — V tanhf^^-1 < 
ICa&E' [ aj l a e JJ 


(4.55) 


2C ar E r E e [ [aj v a « JJ 

The final expression for G in Eq. (4.55) has thermal loading, mechanical 
loading, and t-harmnmftffhani cal interaction contributions. The value of the 
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strain energy release rate is independent of load path. That is, for a given 
load increment, G depends only on the final load state and not on the order in 
which o. and AT are applied. See Appendix A for a proof of the path 
independence of G. 

Substituting Eq. (4.55) into Eq. (4.25) and placing Eq. (4.16) into the 
result, the cracking criterion becomes 


^ E [a.<7. - a,E r (a, - a,)ATf 2tanhf^j-tanhf^ ZG h (4.56) 




2£a r a t E r E, 


Eq. (4.56) must be solved numerically or graphically to find h for a given AT 
and <t.. See Appendix B for a more detailed presentation of both the shear lag 
solution and the derivation of energy expressions. 


4.1.3 Degradation of Laminate Properties 

To include the effects of microcracks in one layer on crack formation in 
other layers, expressions for the reduced laminate properties as functions of 
crack density will be found. Cracks are assumed to reduce the properties of a 
ply group as a function of crack density. Using CLPT, the effective properties 
of the laminate are then determined from the degraded ply properties. Laws 
and Dvorak [24] derived the loss of longitudinal stiffness in a cracked 
laminate. The average strain of the segment between cracks in the 
uncracked portion of a mechanically loaded laminate can be shown to be 




• 1 + 

\ 


a*E 

'tsAzs- 


E 0 1 2C ha r E r 



(4.57) 


This expression can be valid for any two cracks 2h apart. Substituting the 
expression for average crack density, 


P = 


1 

2 h 


(4.58) 
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and rearranging gives the effective stress-strain relation for the cracked 
laminate. 


(4.59) 

where 


EM- 


1 + 


E. 


^tanh 
M 1 



(4.60) 


is the new laminate stiffness as a function of crack density. 

McManus et al. [56] went a step further to derive reduction of all 
laminate properties due to cracking. Considering the reduction in stiffness to 
be caused entirely by a reduction of the effective stiffness of the cracking ply 
group, they define a knockdown factor, k , due to the microcracks 


From Eq. (4.15), 


E e (p)~ *Ec 


(4.61) 


E o x ^E £ + a s E^ 

(4.62) 

Substituting Eq. (4.61) into the equation formed by setting E 0 =E„(p) and 
E c =E e (p) in Eq. (4.62), 


E o (p) , Z& l Bs& s. (4.63) 

Placing Eqs. (4.62) and (4.63) into Eq. (4.60) and solving for k. 


k « a r E r 


1- 



+o <r E e ^tanh 



(4.64) 


This knockdown factor is used to calculate degraded laminate properties due 
to matrix cracks. Details of this method are given in Section 4.2.1. 


I 
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4.2 IMPLEMENTATION 

The implementation of the analysis follows the work of Park [57, 61], 
modified to include thermomechanical loading. The basic formulations for 
crack appearance and reduced laminate properties, Eqs. (4.56) and (4.64), 
respectively, are incorporated into a progressive damage model. This model 
is used to predict crack density and degraded ply and laminate properties at 
incremental loadings of an arbitrary thermo-mechanical load history. It is 
general enough to include any laminate geometry and includes modeling of 
various secondary effects. 

4.2.1 Damage Progression Algorithm 

Assuming an existing crack spacing 2 h, new cracks will form midway 
between existing cracks when the basic energy criterion in Eq. (4.56) is 
satisfied, resulting in a new crack spacing h. An existing crack spacing just 
under 2 h will not satisfy the criteria, and no new cracks will form. In 
practice, the crack spacing is not uniform, and the true crack spacing will fall 
somewhere between these extremes. Hence, the average crack density, p, is: 



where h satisfies Eq. (4.56). For consistency, the present work will use only 
the minimum crack density, where all cracks are separated by 2 h. This 
convention been shown to work well in previous studies [56, 57, 61]. The 
choice is essentially arbitrary, and it will affect the values of G/ e and £ back- 
calculated from the experimental data in Section 4.4. 

The analysis calculates crack density for every ply of any general 
laminate and the resulting degraded laminate properties. The algorithm 
includes effects such as material softening and temperature dependent 
material properties. "Material softening” refers to the fact that cracking in 
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Figure 4.5 Flow char t of algorithm used to implement analysis. 
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one layer will affect the cracking behavior of the other plies in the laminate 
by c h a n g in g the effective properties of the laminate. The reduced properties 
at each temperature or stress increment reflect the damage incurred in all 
plies at all previous loads. Material properties are also known to be 
dependent on temperature. The analysis includes temperature-dependent 
material properties, linearly interpolating between temperatures at which 
data is available. 

Figure 4.5 shows a flow chart describing the basic algorithm. The 
analysis starts at the stress-free temperature and zero mechanical stress and 
increments to the user-specified initial conditions, accounting for any thermal 
residual stresses incurred during manufacture. It then increments through 
the temperatures and stresses in the user-defined load profile. At each new 
load increment, the material properties are obtained from the temperature- 
dependent material property data. A single ply group, referred to as ply c, is 
designated as the cracking layer. The crack density in ply c at this load 
increment is then calculated. 

The analysis is carried out in the local coordinate system of the 
cracking layer, defined in Section 4.1.1. Properties of the rest of the laminate 
are smeared together, whereby existing cracks in the other layers are 
included. Eq. (4.56) is solved for h using a bisection iteration numerical 
method, and crack density p follows from Eq. (4.65). This procedure is 
repeated for every ply group in the laminate at this load increment. Finally, 
knockdown factors for each ply group and the effective laminate properties 
are updated using Eq. (4.64) and CLPT to reflect the total damage of the 
laminate at this load. These steps are iterated to calculate crack density and 
laminate properties through the entire load history. 
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4JL2 Derivation of Effective laminate and Ply Properties 
Laminate theory (CLPT) [63] is used to derive the stiffness constants 
used in the above analysis. The equivalent stiffnesses, E ol E r and E e are 
necessary to solve for crack spacing in Eq. (4.56). First, appropriate 
temperature-dependent material properties are obtained. Each ply i has 
material properties Eu (longitudinal stiffness), E a (transverse stiffness), v< 
(major Poisson's ratio), G< (shear stiffness) ^(longitudinal CTE), and aa 
(transverse CTE). Ply i has thickness t ( . The fibers of each ply are aligned at 
an angle 0 t to the x axis of the global coordinate system. The cracking ply 
group, ply c, is treated as a single layer with orientation de. The crack 
formation analysis is carried out in the local coordinate system x’y'z' of the 
cracking ply group, defined in Section 4.1.1. In this coordinate system, the 
ply angles are defined: 


0/= d i + 9Q°-9 e (4.66) 

The necessary laminate properties for computing crack density and property 
degradation can be calculated using familiar CLPT relations [61]. The 
laminate stiffness in the x'y’z' system is 

A = £Qa (4.67) 

1*1 

The rotated reduced ply stiffnesses in the x'y'z' system are 


Q < =T7 1 Q,T7 t 
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(4.68) 


(4.69) 
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T,= 


coa 2 9' sin 2 9[ 2 sin 6' cos S' 

sin 2 9’ 1 cos 2 -2 sin 9' cos 9' 

- sin 9' cos sin 0/ cos 91 cos 2 0/ - sin 2 9' 


(4.70) 


and the reduced ply stiffnesses are 


q -ElL Q 

W U (i) - Di K i D{ 

Qa«) * ^ ^U> = Ki l} 

" U 


(4.71a-e) 


Here, Ki is the knockdown factor for ply i defined by Eq. (4.64) which accounts 
for the effects of pre-existing cracks. It has the value 1 until ply i begins to 
crack. Note that all matrix-dominated properties are assumed to be reduced 
by the same knockdown factor. The CTEs of each ply are 


a, =< 


<*u 

<** 

0 


(4.72) 


In the xy'z system, the ply CTEs are 


o, = Tfa (4.73) 

The laminate constants required in Eqs. (4.56) and (4.64) are now 
calculated. The total laminate thickness, a., is 

«. = £<, (4.74) 

i-1 

The necessary constants for the cracking layer are 

£, = jf- «, = «. a. = t, (4.75) 

and the smeared properties of the rest of the lamin ate are 

E r = 4nzEs&. a r = a[ a r = a 0 -a c (4.76) 

Or 
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where a[ is the first element of the vector of CTEs which are calculated from 


a r = [A - Q c f e ] | XQi^i - Q c*ct e 


(4.77) 


Eq. (4.56) is now used to predict the crack density in the cracking ply group. 
The knockdown factor for this group is also recalculated using Eq. (4.64). 

After all the ply groups have been analyzed in this manner, the 
degraded effective laminate properties in the global x direction are calculated 
from 


E*. 1 A!"'= A" 1 a* r » a'. (4.78) 

where a* is the first element of 

o' - A-‘ £ 3,04 (4.79) 

iml 

» 

and A, and 5 are calculated from Eqs. (4.67) to (4.73) with $ t = d t . Note 
that bending is not incorporated into this method, so it is strictly valid for 
symmetric laminates only. All of the above steps are repeated for each load 
increment until completion of the entire load history. Each increment 
incorporates the knocked down properties of all the plies from the previous 
increment and temperature-dependent properties for the conditions of the 
current increment. 

4.241 Computer Code 

The computer code CRACKOMATIC II was modified to implement the 
present methodology. Given material properties, laminate geometry, and 
user-defined thermomechanical load history, it predicts crack density and 
corresponding degraded laminate properties. Also included are user options 
to incorporate or omit material softening effects and temperature-dependent 
material properties. The output of the code is a table with columns listing 
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temperature, applied stress, corresponding crack density in any selected 
plies, and effective laminate properties. These can be used to generate plots 
of cracking and changing laminate properties as the laminate is exposed to a 
given load profile. 

The program output reports crack densities as they could be observed 
on the edge of a specimen. The crack density calculated from Eqs. (4.56) and 
(4.65) is expressed in the local x'y’z ' coordinate system. The calculated 
results are transformed to the global xyr- coordinate system, accounting for 
the geometric effect illustrated in Figure 4.6. The reported crack density, p e , 
is the calculated density, p*,, multiplied by a geometric factor 

Pe = Poo sin(0,) (4 >80 ) 

The manual for CRACKOMATIC II is in Appendix C. The document 
describes the program in detail and shows sample sessions. The source code 
is available by request from the author or the TELAC laboratory. 

4.3 LAMINATE SELECTION 

Preliminary analysis was necessary to select the three laminate 
geometries to be used in the experimental investigation. Several criteria 
were used to determine the layups. First, in order to capture both the onset 
and accumulation of cracking, microcracking in the laminates had to initiate 
within the available temperature range of the environmental chamber. If 
cracking initiated at one of the lowest temperature increments or not at all, 
then the test would yield little useful data. Furthermore, the laminates had 
to be initially uncracked at room temperature following cure in order to 
capture the onset of cracking. Using the existing progressive thermal 
analysis of Park [571, it was found that all plies of a [O/454/9O4/-454], l am i n ate 
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Crack density measured on this side 



Figure 4.6 


Higher crack density observed in 90° layer than in 45® layer 
during edge inspection, though crack spacing same in both 
[57]. 
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would indeed begin cracking by -184°C, the lowest temperature possible in 
the chamber, and would be uncracked at room temperature. Table 4.1 lists 
the material properties used in these preliminary calculations. The values of 
G /C and £ were estimated from those of T300/934 graphite/epoxy, found in a 
previous study [56]. 

Based on these preliminary thermal calculations, shown in Table 4.2, 
two additional layups were chosen. A [0 a /45 a /90 a /-45 2 ] a laminate would be 
used to investigate the effects of layer thickness by comparing its cracking 
behavior with that of the [04/45^/90 4/-45J, laminate. Second, a [0a/60a/-60J. 
layup was chosen to examine ply orientation effects by comparison with the 
[Oa/45a/90a/-452], laminate. The progressive thermal analysis indicated that 
the onset of cracking in both of these laminates should occur near the lower 
temperature limit of the environmental chamber, limiting the number of 
temperature increments at which useful crack density data could be collected. 
However, Park's work showed that thin ply groups tend to initiate cracking at 
a higher temperature than predicted by that analysis [57], and it was 
suspected that these two laminates might demonstrate this "thin ply” 
behavior. 

Further preliminary calculations were necessary to check whether the 
three laminates chosen would yield useful data under mechanical loading. 
These are summarized in Table 4.3. In order to maximize the range of 
loading under which microcrack data can be collected, crack initiation could 
not occur just before or simultaneously with laminate failure. CLPT was 
used in conjunction with a maximum stress failure criterion to predict first 
ply failure. A thermal residual stress, assuming a stress-free temperature 
equal to the cure temperature (177°C), was included. 
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Table 4.1 AS4/3501-6 Material Properties Used in Preliminary 

Calculations 




























Table 4.2 


Summary of Preliminary Calculations Using Progressive 
Thermal Analysis from [57] 


Laminate 

Crack Initiation 
Temperature 
(°C) 


-58 


-156 

■m 

-124 
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Table 4,3 Summary of Preliminary Calculations for Crack 

Initiation and Final Failure under Mechanical Loading | 

I 

I 

















Final failure of the laminate was calculated using the ply discount 
method [63]. In this method, CLPT was used with a maximum stress failure 
criterion to predict first ply failure. If the axial failure stress (X c or X t ) was 
exceeded, B h E„ v, and G* for the failed ply were set to zero, corresponding to 
fiber failure. If, on the other hand, the transverse ( Y e or Y,) or shear (S) 
failure stress was exceeded, E„ v, and G u for that ply were set to zero, 
corresponding to matrix failure. The stress was subsequently increased until 
the next ply failed, and the properties of that ply were knocked down as after 
first ply failure. This procedure was repeated until fiber failure in the 0° 
plies is predicted, which was presumed to indicate final failure. The 
predicted behavior of all three laminates, including both first ply failure and 
final failure, was deemed acceptable. 

Next, the preceding analysis was repeated using an in situ transverse 
tensile strength Yp [65], 

Y* = r,(l + A£$®>) (4.81) 

where Ad is the minimum ply angle change between the ply under 
consideration and its neighboring plies, N is the number of consecutive plies 
in that layer, Y* is the transverse tensile strength of a unidirectional ply, and 
A and B are constants of the strength distribution, respectively, for T300/934. 
This empirical formula was originally derived for T300/934, for which A and 
B are usually taken to be 1.3 and 0.7, respectively. It is applicable to 
AS4/3501-6 as well, since the cracking behaviors of the two material systems 
are similar. The predicted first ply failure and final failure loads using this 
in situ transverse tensile strength were deemed acceptable for all three 
laminates. 
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4.4 DATA FIT 

As discussed in Section 2.4, the critical strain energy release rate is 
impossible to measure directly. Many researchers use the interlaminar 
fracture toughness as a value for G Ie . However, as this property is measured 
in double cantilever beam tests, it essentially characterizes delamination 
rather than microcracking. Moreover, measured values of the interlaminar 
fracture toughness vary considerably, as shown in Figure 2.1. The shear lag 
parameter is a semi-empirical aspect of the shear lag model which likewise 
cannot be measured directly. This parameter has not been standardized, 
rendering values used in different studies difficult to compare. 

Consequently, the critical strain energy release rate and shear lag 
parameter will be determined from experimental data using the following 
data fit. First, crack density is predicted as a function of applied load or 
progressive cooling using a particular G h and £. The predicted and observed 
crack densities are normalized by the highest density recorded for that 
laminate, and the applied load or change in temperature is normalized by the 
largest load or AT to which the laminate was exposed experimentally. Next, 
a single error value is calculated for each experimental data point using 


Error j = min 


^TV.-AT.,]' +[p„-p. y f ] 


(4.82) 


for thermal loading. AT«* and are the normalized change in temperature 
and predicted crack density, respectively, at the k* load increment of the 
analysis. and p tJ are the normalized change in temperature and 

observed crack density, respectively, at the/* data point. The value used for 
Error j is the minimum value calculated using any h. Similarly, under 
mechanical loading, 
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Error, = min 




(4.83) 


where Po* and P tJ are the normalized applied load at the k* load increment 
of the analysis and j* data point, respectively. Conceptually, this method is 
equivalent to plotting the normalized experimental data and analytical 
predictions on the same graph and finding the error for each data point from 
the absolute distance between that data point and the nearest point on the 
analysis curve. 

The mean square error (MSE) is then determined for every layer in 
each laminate using [66] 


1 JL 

MSE = — V Error*, (4.84) 

where m is the total number of data points for the layer. This is repeated for 
both thermal and mechanical loading, for each layer in every laminate. A 
simple computer code was written to carry out the calculations. Finally, all 
the MSEs are summed, giving a single total MSE value for that G /e and £. 
The procedure is repeated for different values of the shear lag parameter and 
critical strain energy release rate until the single set of values which 
minimizes the total MSE is found. Those values can then be used in Eq. 
(4.56) to analyze any laminate, provided that the material system does not 
change. 

This data fit has several important advantages over the other methods 
used to determine Gu and £ described in Section 2.4. A large amount of 
experimental data may be reduced to a single set of G lc and £ values. The 
data fit can be used to check for laminate dependence of these parameters by 
finding the variation in values fit for the individual laminates. Finally, if 
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found to be laminate independent, G te and (, only need to be measured once 
for a given material system. 


I 
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CHAPTERS 


EXPERIMENTAL PROCEDURES 


An experimental investigation was carried out to provide greater 
qualitative understanding of the microcracking problem, determine input 
parameters for the analysis, and provide verification of the analytical model. 
In this chapter the specimen matrix used in the investigation, and the 
rationale behind its design, are presented. Manufacturing, pre-test 
preparation, and instrumentation of the specimens are then described. 
Finally, specimen conditioning and data collection procedures are presented. 

5.1 TEST MATRICES 

Small, rectangular graphite/epoxy specimens, referred to as thermal 
specimens, and larger tensile coupons were manufactured using three 
different layups. Thermal specimens were cooled to progressively lower 
temperatures. After reaching each target temperature, specimens were 
returned to room temperature. Microcrack density data was then collected on 
the specimen edges and in their interiors by optical microscopy. Tensile 
coupons were loaded monotonically to progressively higher loads. Edges were 
inspected for cracks after unloading the coupons from each load level. Stress- 
strain data was also collected. 

AS4/3501-6 graphite/epoxy was chosen as the material system for this 
investigation. It is a well understood material, and an extensive database of 
its properties has been established at TELAC. Three different laminates 
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were used: [(V454/9<y-45J„ [0 2 /45 i /90 a /-45 2 ]., and [0a/60a/-60J,. The layups 
were selected based upon the procedure described in Section 4.3. Two 
different sizes of thermal specimens were chosen. The first, 
7.62 cm x 1.27 cm, was identical to that used by Park [57], so that data 
collected might be compared with the results of that study. To investigate 
width effects, a second size was chosen, 7.62 cm x 2.54 cm. The width of the 
tenailft coupons was selected to be the same as these wider thermal specimens 
to allow thermal and mechanical crack density data to be compared. Their 
length was subsequently chosen to be 25.4 cm, based on ASTM D3039. 

A matrix of all of the specimens manufactured is shown in Table 5.1. 
Table 5.2 is a test matrix, which shows the conditioning and examination 
performed on the specimens from a single layup. Specimens made from all 
three layups were tested in the same way. 

5.2 TEST SPECIMEN MANUFACTURE AND PREPARATION 

All laminates were manufactured at TELAC using standard 
procedures of the laboratory [67]. 35.6 cm x 30.5 cm panels were fabricated 
naing standard hand layup techniques. A total of six panels, two of each 
layup, were made. The material for the layup was furnished as 30.5 cm wide, 
continuous AS4/3501-6 unidirectional prepreg tape. Nominal ply thickness 
was 0.134 mm. The panels were autoclave cured according to the 
manufacturer's cure cycle. Following cure, the panels were postcured at 
177°C for eight hours. 

Fewer than twelve hours after postcure, the panels were cut into 
individual thermal and mechanical specimens as per the cutting plan in 
Figure 5.1. The same cutting plan was used for all panels. The cutting plan 
was randomized so that the effects of manufacturing variations, such as 
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Table 5.1 Specimen Matrix 


I 

I 

I 

I 


Layup 

■E53S353Sii HHS3HH BESS^M 


Thermal/Edge 

7.62 x 2.54 

5 



7.62x1.27 

5 

[O 4 / 454 / 9 O 4 M 5 J. 

Thermal/Interior 

7.62 x 2.54 

5 



7.62 x 1.27 

5 


Mechanical/Edge 

25.4 x 2.54 

15 


Thermal/Edge 

7.62 x 2.54 

5 



7.62x1.27 

5 

[OaMSa/SOaMSJ. 

Thermal/Interior 

7.62x2.54 

5 



7.62x1.27 

5 


Mechanical/Edge 

25.4 x 2.54 

15 


Thermal/Edge 

7.62x2.54 

5 



7.62x1.27 

5 

[Oa/ffy-eoj, 

Thermal/Interior 

7.62 x 2.54 

5 



7.62x1.27 

5 


Mechanical/Edge 

25.4 x 2.54 

15 
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Table 5.2 Test Matrix Repeated for All Laminates 




Cool to -18°C, 
return to RT, cool 
to -73°C, return to 
RT, cool to -129°C 
return to RT, cool 
to -184°C return 
toRT 


Tensile loading to 
2.22 kN, 4.44 kN, 
6.66 kN, etc. until 
failure; return to 
no load after each 


Data Collected 


Measure cracks on 
edge at RT before 
conditioning, and at 
RT after each 
target temperature 

5 

Measure cracks in 
interior 

3 

Measure cracks in 
interior 

3 

Measure cracks in 
interior 

3 

Measure cracks in 
interior 

3 

Measure cracks in 
interior 

3 

Measure cracks on 
edge before 
conditioning and 
after each loading 

0 





t Room Temperature (21°C) 























spatial variability in laminate thickness and fiber volume fraction, would not 
create a systematic bias in the results. A water cooled diamond blade was 
used for cutting. 

Both long edges of every specimen were wet sanded with 600 grit 
sandpaper and polished with 0.7 micron grit powder on a polishing wheel. 
The reasons for polishing were to standardize crack initiation sites and to 
facilitate optical microscopy of the specimen edges. Immediately after 
polishing, all specimens were dried at 93°C for 14 hours. They were stored in 
airtight containers with desiccant for the remainder of the experimental 
investigation to minimize moisture absorption effects. 

Some additional preparation was necessary for the tensile coupons. A 
total of 3 width and 6 thickness measurements were taken from each coupon. 
This was done to ensure specimen uniformity and to compare measured 
values with nominal ones. A template of the positions at which thickness and 
width measurements were to be made was aligned to the center of the test 
section of each specimen, and the positions were marked out for 
measurements. Figure 5.2 shows a schematic of a tensile coupon with the 
measuring points. The average width of the coupons was 2.53 cm, compared 
to the no minal width of 2.54 cm. Average ply thickness was found to be 
0.135 mm, compared to the no minal value of 0.134 mm. Standard deviations 
were amall, 0.00292 cm for the width measurements and 0.00422 mm for the 
thickness measurements. All subsequent calculations in the present work 
use nominal values for both width and thickness. 

Next, fiberglass end tabs were bonded to the ends of the mechanica l 
coupons to reinforce the gripping section where the tensile loading was 
applied. The and tabs bad a n ominal thickne ss of 5.3 mm. To ensure smooth 
transfer of loads from the end tabs to the specimens and to minimize stress 
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Figure 5.2 Schematic of tensile coupon measurement locations. 
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8 : 


concentrations, the ends of the tabs nearest to the test section were beveled 
30°. The dimensions of the end tabs were 2.54 cm x 5.08 cm, leaving a 
15.2 cm test section. 

The end tabs were cut from 21 -ply, pre- cured sheets of fiberglass with 
a bandsaw. A belt sander was used to bevel one of the edges of each tab. 
Prior to bonding, both the tabs and the coupons were cleaned with 
cheesecloth. Following a standard TELAC procedure [67], the tabs were 
placed on the coupons at room temperature with a layer of adhesive film 
between the tab and specimen. A total of four tabs were placed on each 
specimen. The adhesive used for bonding was Cyanamid 123-2 Adhesive 
Film. The end tabs were bonded to the coupons using a secondary bond cure 
in an autoclave, carried out at 1.45 mPa gage vacuum and 107°C for two 
hours [67]. 

Tensile coupons were instrumented with EA-06-125AD-120 strain 
gages manufactured by Micro Measurements. The gages were 5 mm x 10 mm 
with a gage factor of 2.080 and an accuracy of ±0.5%. One gage was mounted 
longitudinally on each coupon at the center of the test section, as shown in 
Figure 5.3. After cleaning the specimens thoroughly with cheesecloth, the 
gages were bonded to them using a catalyst and M-Bond 200 adhesive 
supplied by Micro Measurements. Electrical leads, 91 cm in length, were 
soldered to the gages. 

541 TESTING PROCEDURES 

5.3.1 Thermal Conditioning 

Thermal loading specimens were progressively cooled in a thermal 
environment chamber. The chamber used electrical resistance rods for 
heating and liquid nitrogen injection for cooling, allowing a temperature 
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Top View 


Side View 


25.4 cm 



Graphite/Epoxy 


Strain Gage 


Film Adhesive 


Glass/Epoxy 
End Tab 



2.54 cm 

Note: Not to Scale 


Figure 5.3 Tensile coupon configuration. 
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range of -184°C to 427°C. Specimens were placed on stainless steel racks in 
the 30.5 cm x 10.2 cm x 10.2 cm test volume within the chamber. The 
specimens were shielded from direct heat radiation and contact with liquid 
nitrogen and were directly heated and cooled by fan-circulated air only. The 
temperature of the chamber was precisely controlled with an Omega 
temperature controller. The microprocessor-based controller could be 
programmed to a desired thermal profile comprising a series of linear 
segments. A single J-type thermocouple provided feedback to the controller. 
Prior to testing, over 100 tuning runs were done to determine the optimum 
controller tuning settings and feedback thermocouple location. 

Five additional K-type thermocouples were mounted within the 
chamber. One of these was epoxied to a graphite/epoxy specimen to simulate 
the thermal response of an actual test specimen. Thermocouple outputs were 
fed into a MacADIOS breakout box, which was in turn interfaced with an 
Apple Macintosh IIx through a MacADIOS analog/digital converter card 
supplied by GW Instruments. Thermocouple readings were recorded and 
displayed in a virtual strip chart by the LabVTEW®2 data acquisition 
program. The thermocouples were used in the tuning runs to determine the 
controller settings and specimen rack locations which provided stable control 
and minimized the effects of temperature gradients. During thermal testing 
the thermocouples were used to record the actual temperature profile. 

Thermal specimens were first inspected for both edge and interior 
microcracks at room temperature using the procedure described in Section 
5.4. They were then cooled to progressively lower temperatures down to 
-184°C. The target temperatures for all thermal specimens were -18°C, 
-73°C, -129°C, and -184°C. Since multiple runs were necessary to cool all the 
specimens, the individual specimens included in each run as well as their 
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distribution on the specimen racks were always randomized to prevent these 
factors from causing a consistent bias in the results. The randomization 
procedure involved gently shaking the bottle in which the specimens were 
stored and simply pulling them out randomly. Cooling and heating were 
carried out at 13.9°C per minute, slow enough to minimize both rate effects 
and thermal gradients. Figure 5.4 shows the temperature profile of the 
thermal environment chamber during a test run to -73°C, as measured by the 
thermocouple epoxied to a specimen. 

The specimens designated for edge crack data, referred to as "edge" 
specimens, were cooled to the target temperature, soaked at that temperature 
for 5 minutes, returned to room temperature, and soaked at room 
temperature for 5 minutes. Crack density data was subsequently collected on 
their edges. They were then returned to the environment chamber and cooled 
to the next target temperature. This procedure was repeated until all edge 
specimens were conditioned to every target temperature. 

The specimens used for internal crack counting, referred to as 
"interior" specimens, were separated into five different groups for each of the 
five target temperatures. Each interior specimen was conditioned to every 
target temperature down to and including its designated target temperature, 
after which it was destroyed by the inspection procedure. The specimens 
were cooled to a target temperature, soaked at that temperature for 5 
minutes, returned to room temperature, and soaked at room temperature for 
5 minutes. The specimens designated for that target temperature were 
inspected for internal crack density, rendering them useless for further 
testing. The remainder of the specimens were then conditioned to the next 
target temperature, and the specimens designated for that temperature were 
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Temperature (°C) 
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Figure 5.4 Temperature profile of thermal environment chamber during 
test run to -73®C. 
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inspected. This procedure was repeated until all target temperatures were 
reached. 

No more than one day was allowed to pass between conditioning the 
specimens to a target temperature and inspecting them for microcracks. 

5.3*2 Mechanical Conditioning 

The mechanical specimens were first inspected for edge microcracks at 
room temperature using the procedure described in Section 5.4. Tensile 
coupons were then loaded monotonically under load control to progressively 
higher target loads. The coupons were loaded to a target load, unloaded, 
removed from the test machine, and inspected for cracks on their edges. They 
were then placed back in the testing maching, and loaded to the next target 
load, which was 2.22 kN higher the previous maximum load, unloaded, and 
inspected again. This procedure was repeated until other damage modes, 
particularly delamination, were observed. The coupon was then loaded to 
failure, and no further microcrack data was collected. 

All tensile tests were performed on an MTS 810 test machine equipped 
with hydraulic grips. To mount a tensile coupon in the machine, one end was 
first placed in the upper hydraulic grips. Once alignment of the coupon was 
achieved by means of a machinist's square, gripping pressure was applied. 
Next the upper crosshead was lowered until the bottom of the specimen was 
between the lower hydraulic grips. Before gripping the lower jaws, the strain 
gage was balanced and calibrated. Loading and unloading rates were held 
constant at 89 N/s. The strain gage output was fed through Vishay 
Instruments strain gage amplifiers/balances, into the MacADIOS breakout 
box and analog/digital converter card, and recorded by LabVTEW*2. Load 
and stroke information from the test machine was also recorded. 
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5.4 MICROCRACK INSPECTION 

5.4.1 Edge Examination 

Microcracking damage was quantified by optical examination of the 
polished free edge of a specimen. Cracks were counted in a 2.54 cm section in 
the center of one edge of each specimen. The boundaries of the section were 
precisely marked with a white paint marker. A tally counter was used in 
crack counting to minimize error. Each ply group was inspected individually, 
and only microcracks which spanned at least half the thickness of the ply 
group were counted as cracks. Shorter cracks or fiber/matrix debonding were 
often visible but never counted. The crack density recorded was that 
observed on the edge, and does not account for the geometric effect shown in 
Figure 4.6. 

Specimens were randomly assigned for inspection to two different 
people. Each used an identical optical microscope at a magnification of lOOx. 
Cross checks were done by having both persons inspect the same specimen. 
Generally, both observed and recorded the same results. 

5.42 Interior Crack Counts 

Interior specimens were used to determine the crack density within the 
volume of the laminate. These specimens were cooled in the thermal 
environment chamber with the edge specimens. Then their edges were 
inspected using the procedure described in Section 5.4.1. Next, (1) material 
was sanded away from the inspection edge until reaching a desired depth, 
(2) the edge was polished and the 2.54 cm center section marked, and (3) the 
edge was inspected once again. Steps 1-3 were repeated at a total of four 
incremental depths until approximately half of the laminate volume was 
inspected, as shown in Figure 5.5. The 1.27 cm-wide specimens were 
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Section in which 
Cracks Counted 





13.4 mm 


Figure 5.5 Schematic of procedure for interior crack counts. 
7.62 cm x 2.54 cm thermal specimen shown. 






examined at depths of 0 mm, 2 mm, 4 mm, and 6 mm. The 2.54 cm-wide 
specimens were examined at depths of 0 mm , 4 mm, 8 mm, and 12 mm. 

In Step 1, the edge was sanded to within approximately 10% of the 
next depth increment using dry 180 grit sandpaper. Care was taken not to 
allow the specimens to become hot to the touch during sanding. The 
remainder of the material was removed with dry 600 grit sandpaper. The 
width of the specimen was frequently measured with a vernier caliper to 
ensure even, accurate sanding. Following sanding, the edge was polished 
with 0.7 micron grit powder on a polishing wheel. This sanding procedure 
was shown by Park [57] to cause minimal additional damage to the 
specimens. Finally, a 2.54 cm section was marked in the center of the edge, 
and cracks were inspected as in Section 5.4.1. 
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CHAPTER 6 


RESULTS AND DISCUSSION 


Analytical and experimental results are presented in this chapter. 
First, the results of the data fit used to determine G lc and (, are reviewed. 
Mechanical loading results and correlations with analytical predictions are 
then presented. Thermal loading results and correlations follow. A 
parametric study is presented to examine the sensitivity of the analysis to 
different variables. Finally, the results and correlations are discussed. The 
effects of laminate geometry, specimen width, and layer thickness are 
examined. The validity of the analytical model and its assumptions is also 
addressed. 

6.1 DATA FIT 

The results of the data fit are presented in Table 6.1. The values of the 
critical strain energy release rate and shear lag parameter which minimize 
the overall mean square error for all layers are G Ie - 141 J/m? and £ = 1.0. 
Both thermal and machnnirfll loading results were used in this data fit. It is 
important to note that all subsequent analyses in the present study will use 
these values. 

Also shown are tlie Gj e and £ values which minimize the mean square 
error for each of the three layups individually. These were determined by 
fitting the analysis to the mechanical cracking data for each laminate. 
Values fit to the thermal data of the laminates are also 
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Table 6.1 Data Fit Results 


Laminate 

Gfc (J/m a ) 

C 

AllData t 

141 

1.0 


141 

1.1 


141 

1.2 


120 

1.2 

Thermal^ 

156 

0.6 


t Fit to all data available; subsequently used in all analyses, 
tt Fit to data from [O 4 / 464 / 9 O 4 /- 454 ]. thermal specimens 
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shown. Thermal results of the [0 a /45 3 /90 a /-45 2 ], and [(V6<V-60 a ], laminates 
were not included in the data fit, as very little cracking was observed in those 
specimens under progressive cooling. The material properties used in the 
data fit and every analysis thereafter are shown in Table 6.2. These values 
were assumed to be independent of temperature. 

The G lc and £ values fit to individual layups and loading conditions are, 
in general, very close to those fit to all measurements. G Ie fit to data from 
mechanically loaded [0a/60a/-60 2 ]. specimens is approximately 15% lower, and 
Gic fit to data from thermally loaded [(V454/9(V-45J, specimens is 10% 
higher, than that fit to all data. The only severe disagreement is the £ fit to 
data from thermally loaded specimens. 

6.2 MECHANICAL LOADING RESULTS 

6.2.1 Crack Density Data and Analytical Predictions 

In this subsection, selected analytical and experimental crack density 
correlations for the tensile coupons are presented. Average measured crack 
densities are plotted with error bars representing one standard deviation. All 
such results not presented here may be found in Appendix D. 

Microcracks observed on the coupon edges generally appeared 
continuous through the thickness of the ply group and were nearly 
perpendicular to the coupon sides. Some cracks were observed which formed 
less than halfway through the thickness of the ply group. These not were not 
included in the reported crack density because they did not meet the counting 
criterion defined in Section 5.4.1. At higher load increments, cracks were 
often observed to branch off of existing cracks. These were only counted if the 
"branch” formed across at least half of the layer thickness. Some of these 
formed at very shallow angles to the xy plane. Figure (6.1) shows one such 
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Table &2 AS4/3501-6 Material Properties Used in All Analyses 


E, (GPa) 

142.0 

a,(p£/° C) 

-0.36 

E, (GPa) 

9.81 

C) 

28.8 

V 

0.30 

T a (°C) [68] 

177 

G u (GPa) 

6.0 

■RBI 

0.134 
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Figure 6.1 Photomicrograph of edge of [04/45^/90 4/-45J. tensile coupon at 
100x magnification. Local delamination propagates from 
shallow-angle cracks visible in 90 4 layer. 
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crack. Local delamination propagated from these kinds of cracks at the 
highest load levels. 

Analytical and experimental results for the 45 4 ply group in the 
mechanically loaded [0 4 /45 4 /90 4 /-45J i specimens are shown in Figure 6.2. 
Correlation between experiment and analysis is very good. The layer has a 
s mall number of cracks before any conditioning. The analysis predicts the 
onset of extensive cracking nearly perfectly and generally falls within the 
error bars during crack accumulation. Scatter is relatively low, except in the 
region where cracking initiates. 

The results for the 90 4 ply group in the same laminate are shown in 
Figure 6.3. The crack density rises to a significantly higher value than in the 
45 4 ply group, though both show the same trends. Cracks initiate at a 
slightly lower load than predicted. The analysis also underpredicts crack 
density through the entire range of cracking. A slight increase in the slope of 
the crack data is visible at approximately 30 kN. This load coincides with the 
first observation of microcracks at very shallow angles during microscopic 
inspection. Local delamination was observed to propagate from these cracks 
at higher load increments, as shown in Figure 6.3. The increase in slope also 
coincides approximately with the onset of cracking in the adjacent 45 4 ply 
group. 

Correlation between experiment and analysis is not as good for the -45$ 
layer, shown in Figure 6.4. As predicted, fewer cracks form in this layer than 
in the 90 4 and 45 4 layers. However, crack initiation occurs at a much higher 
load than predicted. After the onset of cracking, the analysis captures the 
cracking trend satisfactorily. A comparison of Figure 6.4 with Figure 6.2 
shows that cracking initiates in the middle -45$ layer at exactly the same load 
as in the adjacent 90 4 layer. Additionally, an increase in the slope of the 
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Figure 6.2 Experimental results and analytical predictions of crack 
density vs. progressive applied load. 40i ply group or 
[(V45/904/-45J, laminate. 
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[0 745/907-45 1 

4 4 4 4 8 



Figure 6.3 Experimental results and analytical predictions of crack 
density vs. progressive applied load. 90 4 ply group or 
[O 4 / 464 / 9 O 4 /- 454 ], laminate. 
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[0/457907-451 

4 4 4 4 8 



Figure 6.4 Experimental results and analytical predictions of crack 
density vs. progressive applied load. -45$ ply group ot 
[O 4 / 454 / 9 O 4 /- 45 J, laminate. 
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crack density data in the -45 8 ply group can be seen at about 27 kN, the load 
at which cracking initiates in the 45 4 ply group. 

The applied stress at which cracking occurs in the 16-ply 
[(V45 a /90 a /-45 2 ], specimens is double that in the 32-ply specimens. Scatter is 
relatively high in the region of crack initiation. Once cracking initiates, the 
predicted and measured crack densities reach much higher values than in the 
32-ply [(V45/904/-45J. specimens. Anomalous behavior in the 
[0 a /45 a /90 a /-45 a ]« tensile coupons shows nearly identical trends with the 
[O 4 / 454 / 9 O 4/-45 J, laminate data. This resemblance is readily seen in the 45* 
and 90 4 layer data in Figure 6.5. The -45 4 ply group shows the same middle- 
layer behavior discussed earlier, shown in Figure 6.6. Crack initiation in this 
ply group coincides with initiation in the adjacent 90 4 layer. A slight increase 
in the observed crack density slope occurs simultaneously with the onset of 
cracking in the 45 4 layer. 

Analytical and experimental correlation for the 60 2 layer in the 
[0 a /60 a /-60 2 ], tensile coupons is shown in Figure 6.7. The analysis captures 
the onset of cracking but slightly underpredicts the accumulation of cracks. 
The applied stress at which cracking occurs is comparable to that in the 16- 
ply [0 a /45 2 /90 a /-45 a ]. laminates. The middle-layer crack suppression discussed 
earlier is readily apparent in the -60 4 layer data, shown in Figure 6.8. 
Cracking in this ply group, though predicted at a much lower load, initiates 
simultaneously with the adjacent layer. After the onset of cracking, the 
analysis captures the accumulation of cracking nearly perfectly. 

6.2J2 Stress-Strain Data 

In the present investigation, the first observation of extensive 
microcracking during edge inspection will be referred to as first ply failure. 
When the laminates were conditioned at levels above the first-ply-failure 
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Figure 6.5 


Experimental results and analytical predictions of crack 
density vs. progressive applied load. 45 a and 90 a ply groups of 
[0a/45^90a/-45j, laminate. 
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[0 2 /45 2 /90 2 /-45 2 ] a 



Figure 6.6 Experimental results and analytical predictions of crack 
density vs. progressive applied load. -45* ply group oi 
[0^45^90a/-45aL laminate. 
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Figure 6.7 Experimental results and analytical predictions of crack 
density vs. progressive applied load. 60 8 ply group of 
[0»/60y-60s], laminate. 
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[0 2 /60/-60 2 ] 8 



Figure 6.8 Experimental results and analytical predictions of crack 
density vs. progressive applied load. -60 4 ply group of [0a/6(y- 
60J, laminate. 
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load, cracking and popping noises were clearly audible. These acoustic 
emissions presumably coincided with the formation of new microcracks. They 
were only heard when the load exceeded the target load of the previous 
loading. Reloading to the previous maximum load caused no further cracking 
or popping noise. This behavior indicates that unloading and reloading to 
previous target loads had little effect on the internal damage. 

The stress-strain curve of the [0^458/90^-45^, tensile coupons loaded to 
15.6 kN, below its first-ply-failure load, is shown in Figure 6.9. The loading 
and unloading paths appear straight and nearly identical. No cracks were 
observed during edge inspection after unloading. After conditioning this 
coupon to the next target load, 17.8 kN, extensive cracking was observed 
during edge inspection, indicating first ply failure. The stress-strain curve 
for this target load is shown in Figure 6.10. 

An abrupt decrease in slope of this stress-strain curve can be seen at 
300 MPa, which is the peak stress of the previous load increment, shown in 
Figure 6.9. The change in slope corresponds to first ply failure and indicates 
a loss of stiffiiess due to microcracking. It also coincides with the acoustic 
emissions described earlier. Superimposing Figures 6.9 and 6.10 shows that 
the loading region of the stress-strain plot for the 17.8 kN target load has the 
same slope as both the loading and unloading regions of that for the 15.6 kN 
target load. The unloading region of the 17.8 kN stress-strain plot has a 
slightly lower slope than the loading region, indicating a loss of stiffiiess due 
to microcracking. All stress-strain curves for target loads where first ply 
failure was observed may be found in Appendix E. First-ply-failure loads are 
summarized in Table 6.3. 

After first ply failure, subsequent loadings produced cracking and 
popping noises when the load exceeded the target load of the previous 
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Figure 6.9 Stress-strain plot for FA-2M [<y45a/90a/-45 a 
progressively loaded to 15.6 kN. 


8000 
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loading. However, no further anomalies were observed in the stress-strain 
curves during these loadings. All of the stress-strain curves for loadings 
beyond first ply failure were linear to within the resolution of the data. 

6^3 Ultimate Failure Stress 

Ultimate failure stress data for the three laminate geometries is 
presented in Table 6.3. The [0^45^90 jMSa]. specimens failed at an average 
of 592 MPa with a standard deviation of just 2.77 MPa. The [0a/6<y-60 3 ], 
specimens also had very low variability in failure stress. The failure stresses 
of the [0JASJ90J-A5 J, coupons varied considerably, however. These failed at 
an average of 453 MPa with a standard deviation of 40.1 MPa. 

It is also significant to note that the failure stress of the 
[O 4 / 454 / 9 O 4 /- 45 J. coupons was 32% lower than that of the [02/45a/90a/-45J, 
coupons. The only distinction between the two laminates was layer 
thickness, which caused their microcracking behavior to be different. 
Therefore microcracking probably played a significant role in the failure 
mechanisms, dominating the failure of the coupons. 

6.3 THERMAL LOADING DATA 

6.3.1 Interior vs. Edge Data 

As discussed in Section 5.4, some thermal specimens were used to 
collect interior crack data at various depths. This data was subsequently 
averaged over the measurement depths. Other thermal specimens were used 
for edge crack data only. Very little disparity between the interior and edge 
data is observed. This is in contrast to previous studies [57, 61], where 
interior and edge data were significantly different, and crack density varied 
erratically with measurement depth. The interior data in the present work 
shows very little variation of crack density with depth. 
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Table 6.3 First Ply Failure and Ultimate Failure Data 



First-Ply-Fa 

(Ml 

lilure Stress 
PA) 

Ultimate Fa 

(m: 

lilure Stress 
Pa) 

Laminate 

Average 


Average 



138 

19.3 

453 

40.1 


319 

38.5 

592 

2.77 


368 

31.5 

666 

14.2 
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Cracks observed on the specimen edges and in their interiors had the 
same appearance as those on the edges of the tensile coupons, described in 
Section 6.2.1. New cracks were observed to branch from existing cracks 
following conditioning to the lowest target temperatures. Delamination was 
never observed during the inspection of thermal specimens. 

6&2 Crack Density Data and Analytical Correlations 
In this subsection, selected correlations between analytical and 
experimental crack densities are presented. Average measured crack 
densities are plotted with error bars representing standard deviation. All 
thermal loading results may be found in Appendix D. 

Crack data and predictions for the 45 4 ply group in the thermally 
loaded [(> 4 / 454 / 90 / 454 ], interior specimens are shown in Figure 6.11. The 
interior of this layer is initially uncracked at room temperature. The analysis 
appears to successfully capture crack initiation. Agreement at lower 
temperatures as cracks accumulate is also very good. The predicted crack 
density falls within the error bars at all but the lowest temperature 
increment, where the analysis slightly overpredicts the amount of cracking. 
Note that the crack density under thermal loading in this layer is almost the 
same in magnitude as under mechanical loading, shown in Figure 6.1. The 
variability is slightly higher in the thermal specimens, however. 

The results for the 90 4 ply group in the same laminate are shown in 
Figure 6.12. The crack density rises to a slightly higher value than in the 45 4 
ply group, though both show the same trends. The analysis captures crack 
densities at the lowest temperature increments. The laminate is initially 
cracked at room temperature, obscuring the temperature at which cracking 
initiates. It appears that new cracking initiates at the -73.3°C temperature 
increment, which the analysis appears to capture fairly well. Thermal 
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[0/45/90/-45 1 

4 4 4 4 s 



Figure 6.11 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 45 4 ply 
group of [0y45y90i/45 J. interior specimens. 
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Figure 6.12 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 90* ply 
group of[<V45y904/-45J. interior specimens. 
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initiation is more dearly seen in the edge data, shown in Figure 6.13, which 
is very similar to the interior data but shows fewer initial cracks. The 
analysis clearly predicts thermal initiation but slightly overpredicts crack 
density at the lowest temperature increment 

Correlation between analysis and experiment for the -45 8 layer, shown 
in Figure 6.14, shows the same trends as the mechanical loading results. As 
predicted, fewer cracks form in this layer than in the 90 4 and 45* layers, a 
result of the thickness effects discussed earlier. The analysis seems to 
capture crack accumulation well. Further cracking (from the initially cracked 
state) appears to initiate at -73.3°C, the temperature at which extensive 
cracking was first observed in the adjacent 90 4 layer. However, the 
temperature at which cracking initiated in the -45 s layer is difficult to assess, 
due to the limited number of data points available and the presence of cracks 
prior to conditioning. 

Few cracks are observed in the other two laminates under thermal 
conditioning. Analytical and experimental correlations for the 90 3 and 45s 
layers in the [0^45a/90^-452], specimens are shown in Figure 6.15. The 
analysis predicts the onset of cracking in both layers nearly perfectly. Edge 
data and interior data are similar for both layers. Cracking in the -45 4 ply 
group of this laminate, shown in Figure 6.16, begins at a lower temperature 
than predicted. Crack initiation occurs at the same temperature increment 
as in the adjacent 90§ layer. Edge and interior data for the -45 4 layer are 
almost indistinguishable. All crack densities in the [OaMSa/SOaMSal. 
specimens under thermal loading are significantly lower in magnitude than 
under mechanical loading. 

Experimental results and analytical predictions for the 60 a and -60 4 
layers are presented in Figure 6.17. The analysis predicts cracking in both 
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[ 07467907 - 45 .] 

4 4 4 4 8 



Temperature (°C) 


Figure 6.13 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 90* ply 
group of [O/454/9O4/-45J. edge specimens. 
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4 4 4 4 s 



Temperature (°C) 


Figure 6.14 Experimental results and analytical predictions of 
density vs. progressively decreasing temperature. -4 
group of [0/45790 7-45 J, edge specimens. 
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Temperature (°C) 


Figure 6.15 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 90 a and 45 a 
ply groups of [0^45a/90y-45 2 ], interior specimens. 
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[0 2 /45 2 /90 2 /-45 2 ], 



Figure 6.16 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. -45 4 ply 
group of [0^45^90/452], interior specimens. 
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[ 0 2 / 60 2 /- 60 2 ] s 



Temperature (°C) 


Figure 6.17 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. -60* and 60 a 
ply groups of [(V6O3/-6O2], edge specimens. 
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layers within the temperature range of the environmental chamber. 
However, no microcracks are observed in any of the [0a/60a/-60a], thermal 
specimens during progressive cooling. In specimens conditioned to the lowest 
temperature increment, -184°C, very small microcracks are visible at the 
interface between the 60s and -60 4 layers under 200x magnification. These 
cracks span much less than half the layer thickness. Since they would have 
had to span at least half the layer thickness to be counted as microcracks, as 
per the counting criterion given in Section 5.4.1, these cracks are not included 
in the crack density data. 

6.3.3 Width Effects 

Thermal specimens were cut to two different widths, 1.27 mm and 
2.54 mm, to investigate width effects. Specimen width does not appear to 
correlate with any of the observed cracking behaviors. The crack densities 
measured for both the narrow and wide specimens for each target 
temperature appear to have very similar distributions. Moreover, internal 
inspections show that crack density is virtually constant with measurement 
depth in all interior specimens, regardless of width. Thus no width effects 
are observed in the present investigation. 


6.4 PARAMETRIC STUDY 

The sensitivity of the analysis to the input parameters is presented in 
this section. The parameter under examination is varied while all others are 
fixed at the values used in the analytical correlations, given in Tables 6.1 and 
6.2. The effects of the properties are plotted for the 45 4 layer of a 
[O 4 / 454 / 9 O 4 /- 45 J, laminate. 
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The effect of longitudinal ply stiffness E, on predicted crack density 
curves under mechanical loading is shown in Figure 6.18. Increasing 
stiffness has the effect of raising the laminate's resistance to microcracking. 
Crack initiation occurs later, and microcracks accumulate at a lower rate. 
Under mechanical loading, transverse microcracking in a layer is suppressed 
by the constraint of adjacent layers. If the adjacent layers are stiffer, they 
will provide more constraint, explaining the trends in Figure 6.18. Stiffer 
longitudinal plies also carry more of the applied stress. This lowers the stress 
in the transverse plies, and fewer microcracks form. Completely different 
behavior is seen under thermal loading. Here the adjacent layers provide 
constraint, but in this case the constraint causes mechanical stresses to arise 
within the laminate when it is exposed to a change in temperature. It is 
these mechanical stresses which give rise to microcracks. Thus the analysis 
predicts more cracking as E, is increased under thermal loading, as shown in 
Figure 6.19. The analysis is more sensitive to E t under mechanical loading 
than thermal loading. 

If the transverse ply stiffness E t is increased, more of the applied 
mechanical load will be carried by the transverse plies. Under thermal 
loading, increasing the transverse ply stiffness will raise the stresses due to 
constrained thermal deformations. As a result, cracking initiates more 
easily, under both thermal and mechanical loading, when E t is increased. 
This effect is predicted by the analyses shown in Figures 6.20 and 6.21. The 
analysis curve is shifted horizontally with only a very slight change in shape. 
The sensitivity of the analysis to E t under both types of loading appears to be 
about the same. 

Similar trends are seen with critical strain energy release rate G Ie . 
This parameter is a measure of fracture toughness by definition, so one 
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Figure 6.18 Effects of longitudinal ply stiffness E t on predicted cracking 
behavior under mechanical loading. 45* ply group of 
[O 4 / 454 / 9 O 4/-45 J. laminate. 
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[0746/907-46 J 
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Temperature (°C) 


Figure 6.19 Effects of longitudinal ply stiffness E t on predicted cracking 
behavior under thermal loading. 45* ply group of 
[ 07464/90 4/-45J, laminate. 
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[07457907-45 1 
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Figure 6.20 Effects of transverse ply stiffness E t on predicted cracking 
behavior under mechanical loading. 45* ply group of 
[O 4 / 454 / 9 O 4 /- 45 J. laminate. 
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Figure 6.21 Effects of transverse ply stiffness E t on predicted cracking 
behavior under thermal loading. 45* ply group of 
[0 4/45y90 4/-45J. laminate. 
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expects that the laminate's resistance to fracture will rise when G Ie is 
increased. The analysis curve shifts horizontally, and its shape changes only 
slightly. Thus cracking initiates at a higher load or lower temperature, but 
the accumulation of cracks occurs at about the same rate. The effect is the 
same for mechanical and thermal loading, shown in Figures 6.22 and 6.23, 
respectively. Note that the analysis in some cases predicts initial cracking at 
room temperature, due to the 4T=170°C from the stress free temperature. 

Layer thickness has an interesting effect on predicted cracking 
behavior. Predicted crack densities are plotted in Figure 6.24 for the 45, 
layer in a [0^45^90^-45,], laminate as n varies from 1 to 4. There is little 
difference in the n=3 and n* 4 curves. However, decreasing n to 1 increases 
the initiation load significantly. The rate at which cracks accumulate 
increases significantly as n is decreased. The crack density curve shifts 
horizontally, and its shape changes as well. The same trends are seen under 
progressive cooling, shown in Figure 6.25. 

Figures 6.26 and 6.27 show the effect of varying the shear lag 
parameter £. Increasing £ delays the onset of cracking. The slope of the 
crack density curve is increased, so cracks accumulate at a higher rate. This 
effect is similar to the effect of decreasing the layer thickness, and it is 
observed under both thermal and mechanical loading. 

6.5 DISCUSSION 

In thia section, the results and correlations presented in Sections 6.1 
through 6.3 will be discussed. The assumptions used to derive the analysis 
will be critically examined. 
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Figure 6.22 Effects of critical strain energy release rate G fc on predicted 
cracking behavior under mechanical loading. 45* ply group of 
[O 4 / 454 / 9 O 4 /- 45 J, laminate. 
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Figure 6.23 Effects of critical strain energy release rate G fc on predicted 
cracking behavior under thermal loading. 45 4 ply group of 
[0^45/90 4/-45J, laminate. 
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Figure 6.24 Effects of layer thickness on predicted cracking behavior under 

mechanical loading. 46, ply group of f0,/4 5^90^-45 J, 
laminate. 
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Figure 6.25 Effects of layer thickness on predicted cracking behavior under 
thermal loading. 45, ply group of [0^45^90^-45 J, laminate. 
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Figure 6.26 Effects of shear lag parameter } on predicted crad^ belmvior 
under mechanical loading. 45* ply group of 
laminate. 
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Figure 6.27 Effects of shear lag parameter £ on predicted cracking behavior 
under thermal loading. 45 4 ply group of [O 4 / 454 / 9 O 4 /- 45 J. 
laminate. 
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6.5.1 Data Fit 

The fact that the fitted G Ie and C values for the [O4/454/9O4/-45J, and 
[03/45^90 J/- 45 J], laminates are nearly identical is a strong indication that G k 
and £ are independent of layer thickness. The shear lag parameter found for 
the [ 0 a/ 60 ^- 60 j. specimens is the same as that for the [ 0 2 / 45 y' 90 a /- 45 2 ], 
specimens. Thus £ appears to be insensitive to laminate geometry. The fitted 
Gu value for the [(V6O3/-6O2], data shows only a slight disagreement with that 
for the [(V 45 */ 9 <V- 45 J. and WJAZJWJ-A5& data. 

The value of G /e fit to the thermal data from the [(V454/9(V-45J, 
specimens agrees well with the value fit to all data. The shear lag parameter 
fit to the thermal data was significantly lower than that fit to all data. This 
is due to the fact that cracks were observed to accumulate at a slightly lower 
rate with thermal conditioning. This may indicate temperature dependence 
of the shear lag parameter, which is a function of the layer moduli. However, 
as the crack density data used for this fit was limited to a single laminate 
with relatively few data points, the results of the thermal data fit are 
inconclusive. 

The G tc and £ values found by fitting to data from the individual layups 
in the present study show very little variation. The data shown here 
supports the view that the critical strain energy release rate and shear lag 
parameter are material properties, although they may be temperature 
dependent. This justifies the subsequent use of the same G Ie and £ values for 
all analyses. 

6.5.2 Crack Density Results and Correlations 

Nearly all of the [O4/454/9O4/-45J, tensile coupons and thermal 
specimens have cracks in some layers at room temperature prior to any 
thermal or mechanical conditioning. The analysis predicts this laminate to 
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be uncracked at room temperature following cure. A possible explanation for 
this is that the laminates were cooled too quickly at the end of the cure cycle. 
They may also have incurred damage when cut into individual specimens. 
These explanations are supported by the fact that the crack density of these 
specimens remains unchanged for severed load or temperature excursions. If 
crack initiation due to thermomechanical loading had been reached at room 
temperature, one would expect further cracking in these specimens after the 
first loading. Furthermore, the load or temperature at which further 
cracking from this initially damaged state is observed usually coincides with 
predicted initiation. 

Decreasing the thickness of the cracking layer has an effect seen in the 
crack density data of the present study and in many other investigations. 
The layer appears to become more resistant to cracking, since initiation is 
delayed to a lower temperature or higher applied load. Decreasing the layer 
thickness increases the rate at which cracks accumulate after initiation, as 
shown by the parametric study in Section 6.4. 

Very little difference was observed between the crack density data 
collected from the interiors and edges of the thermal specimens. Crack 
density was virtually constant across all measurement depths in the interior 
specimens. Moreover, width effects did not appear to be significant. These 
observations are consistent with one of the most important assumptions of 
the analytical model, that cracks are continuous through the specimen width. 

In every laminate in the present investigation, the onset of cracking in 
the middle ply group coincides with crack initiation in adjacent layers. This 
phenomenon is clearly visible under mechanical loading; thermal results 
were less conclusive. One of the key assumptions of the analytical model is 
that critical starter cracks exist in the cracking layer from which new cracks 
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may nucleate. These defects should be on the order of one ply group 
thickness in extent. In the thicker middle layers, this starter-crack condition 
may not satisfied until the adjacent layer cracks. Then cracks in the adjacent 
layer cause local stress concentrations which act as effective starter cracks. 

There appears to be a correlation between the first observation of 
shallow-angle microcracks in the 90 4 ply group of the tensile 

coupons and the abrupt increase in slope of the cracking data. As shown in 
Figure 6.3, local delamination later propagated from these cracks. Thus the 
appearance of the shallow-angle cracks may have coincided with the onset of 
other damage modes, particularly delamination. The change in slope may 
also be explained by the fact that it coincides with crack initiation in the 
adjacent 45 4 ply group. The microcracks in the 45 4 layer cause local stress 
concentrations at the interface between the ply groups. These may behave as 
initiation sites for both further cracking in the 90 4 layer and delaminations 
between the ply groups. Similar behavior is observed in the [0 2 /45 a /90 a /-45 2 ], 
tensile coupons. The rate of cracking in the -45 4 ply group increases slightly 
at the load increment where cracking initiates in the 45 2 ply group. 

The two behaviors described above are not represented by the 
analytical model. The model assumes that critical starter cracks exist which 
act as crack nudeation sites. If these starter cracks are not present, 
initiation will be suppressed. The model also assumes that cracks in different 
ply groups interact through material softening only, described in Section 
4.1.3, and does not consider local stress concentrations. Thus it does not 
capture some of the observed interactions. However, the analysis is 
conservative in the sense that it always assumes that defects exist. If these 
defects are not present, then crack initiation is suppressed, and the analysis 
gives a conservative prediction of crack density. After the onset of cracking, 
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the starter-crack assumption is satisfied, and the observed cracking behavior 
again agrees with that predicted by the model. This effect can be seen rather 
dramatically in Figure 6.8. 
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CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 


The present work was conducted to investigate microcracking in 
composites. The primary research goal was to develop an analytical 
methodology to predict microcracking in composite laminates under thermal 
and mechanical loading. This goal was met through a combination of 
analytical modeling and experimental investigation. In this chapter, 
conclusions are drawn from the findings of the present work, and directions 
for future research are recommended. 

7.1 CONCLUSIONS 

The discussion presented in Chapter 6 leads to the following 
conclusions: 

1. The primary objective of the present investigation was met. 
Correlation between analysis and experiment was generally very good, 
verifying the analytical methodology. 

2. The same analytical model successfully predicts cracking behavior 
under both progressive cooling and mechanical loading. This indicates 
that the cracking mechanisms are similar for each and may be 
captured by a unifi ed analytical model. Though no combined 
thermomechanical loading was used, whereby progressive cooling and 
mechanical loading occur simultaneously, the mechanical loading 
included thermal residual stresses due to manufacture. 


137 


3. 


The shear lag parameter and critical strain energy release rate did not 
show laminate dependence, verifying the assumption that they are 
material properties. This has important implications for the future use 
of the methodology, as these parameters need to be measured only once 
for a given material system. This may be done relatively simply by 
loading a cross-ply laminate under monotonic tension, collecting crack 
data at progressive loads, and using the analysis to back out the 
parameters. 

4. The laminates in the present work did not demonstrate the "thin ply" 
behavior observed in previous studies [57]. Crack densities were 
similar at specimen edges and throughout specimen interiors, and no 
dependence on specimen width was seen. These results are consistent 
with one of the most important assumptions of the analytical model, 
that cracks are continuous through the specimen width. 

5. Ply groups appeared to interact with one another through local stress 
concentrations caused by cracking in adjacent layers. This behavior is 
not predicted by the analytical model, which is global in nature. 

6. The onset of microcracking was suppressed in middle layers, 
presumably due to a lack of crack nucleation sites. The analysis 
assumes that critical starter cracks always exist, which may not have 
been true in these cases. Cracking initiated instead after adjacent 
layers began to crack, whereby local stress concentrations probably 
acted as nucleation sites. 

7. The analysis is inherently conservative in the sense that it assumes 
that damage is always present. In the case where cracking was 
suppressed in the middle layers, the analysis instead predicted that 
cracking would occur. After the observed initation of cracking, the 
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presence of damage is consistent with the analytical model, and the 
observed behavior nearly matches the predicted behavior. 

7J2 RECOMMENDATIONS 

The present work raises a number of issues which warrant further 

investigation: 

1. Further data collection is necessary, including 

a) more thermal data, over a broader temperature range and with more 
closely spaced target temperatures, to provide better verification of the 
methodology under thermal loading and assess the temperature 
dependence of the material properties; 

b) true thermomechanical loading to verify the thermal and mechanical 
coupling predicted by the analysis and to confirm the path 
independence of the strain energy release rate; and 

c) data for a larger number of laminates to verify the laminate 
independence of G Ic and 

2. Some interesting cracking mechanisms should be explored, such as 

a) whether the analysis has the ability to capture cracking behavior in 
"thin” ply groups (as defined by Park [57]) and 

b) whether the observed suppression of crack initiation in the middle ply 
groups is a general phenomenon. 

3. More complicated global-local models may allow the present 
methodology to predict nearly all of the observed behavior. This 
includes modeling of 

a) local stress concentrations to capture the anomalous interactions 
between adjacent layers, 

b) crack nudeation sites to predict the suppression of crack initiation, and 
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c) interaction between microcracking and delamination to allow the 
prediction of failure dominated by these mechanisms. 

4. The method could be applied to laminates made with materials other 
than AS4/3501-6 by collection of G le and £ values for these materials. 
G/ e and £ appear to be material properties. They could be collected by 
fitting to data from simple tests, such as mechanically loaded cross-ply 
laminates, using the analysis presented here. 
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APPENDIX A 


PATH INDEPENDENCE OF STRAIN ENERGY 

RELEASE RATE 


A brief proof of the path independence of the strain energy release rate, 
G, is presented in this section. From Eq. (4.54), 


G = ip[a 0 <x a - *rEr ( a c - a r ) ATf <3> 


where 


= , q A. _ 

2Ca 0 a r E r E. 


(A.1) 


(A.2) 


1 

N 

- tanh 

r 2 CftY 

L J 


, °C J. 


Taking the partial derivative of Eq. (A.1) in /, the load increment path, 


(A.3) 


SG 

SI 


^ ~ 2a t a r E r {a c + AT 

+ 2a r 2 E?(a e -a r ) 2 Ar^j<I» 


(A.4) 


Three different load paths (A, B, and C) are shown in Figure A.1. The initial 
and final load states for all three are points 1 and 3, respectively. The stress 
and change in temperature at the final load state are o/ and AT/>, respectively. 
The strain energy release rate at point 3 is found by integrating Eq. (A.4) 
along line /. 

Path A is a single continuous line from points 1 to 3 given by 
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<j a = tjAT 

Substituting Eq. (A.5) into Eq. (A.4) and integrating along 7, 

G a = vj JT' 2a 0 a <J a «5<T a - 2a 0 a r E r (a c - a r )j£' 2^S<r a 
+£' 2a* E* (a e -a r ) 2 ^ 8a a j<D 


(A.5) 


(A.6) 


From Eq. (A.6), 


„ ( a r E r (a e -a r )) A 

G A = r o>l a, ^ L j * 


(A.7) 


Path B is piecewise continuous, represented by a line from points 1 to 2 
and one from points 2 to 3. These are given by 


and 


respectively, where 


<T o =0 


AT = AT f 


AT r = — 


(A.8) 


(A.9) 


(A. 10) 


Substituting Eqs. (A.8) to (A. 10) into Eq. (A.4) and integrating along 7, 

G 8 = v^jj; 2o r 2 E r a ( a e - a r f ATSAT d> 

+ 2a?<7.5<7 4 -£' 2a g a r E r ( a e - a r ) -^ 5a a j<I> 


(A. 11) 


which gives 




(A. 12) 
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Path C is piecewise continuous, represented by a line from points 1 to 
2' and one from points 2' to 3. These are given by 

AT = 0 (A. 13) 

and 

<7 a = Of (A. 14) 

Substituting Eqs. (A. 13) and (A. 14) into Eq. (A.4) and integrating along I, 
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Path (A) 


Path (B) 


Path (C) 


3 





Figure A.l Three load paths used in proof. 


APPENDIX B 


DERIVATION OF SHEAR LAG SOLUTION 
AND ENERGY EXPRESSIONS 


B.1 SHEAR LAG SOLUTION 

From equilibrium of the laminate, 


CT a a o = <V*r + <?c a c 

From equilibrium of the cracked layer, 


-<*P e + 


^c + ^f-d*’)o e -2q<k , = 0 


and solving for q. 


v 2 d*' 


From equilibrium of the rest of the laminate, 


and solving for q, 

a, dc r 
a - — c. — c. 

H 2 dx' 

The stress strain relation for the cracking layer is 

a c = E e (e e - ct c A T) where £ c = ^- 
so 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B.6) 
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alT 

E, Ax' * 

The stress strain relation for the rest of the laminate is 


(B.7) 


a r = E r {e r - a r AT) where e r = 


(B.8) 


o, du, at 

t v' ’ 


(B.9) 


q = K(u e -u r ) (B.10) 

K is a s tiffile ss constant relating the displacements of the two layers to the 
shear stress between them. It can be determined from 

K = — (B.ll) 

<V 

where (X is the effective shear modulus of the shear transfer region. Placing 
Eq. (B.ll) into Eqs. (B.3) and (B.5), 


Kiu — u ) = — — 
r) 2 dx’ 

(B.12) 

K(u u ) = - flr dffr 
' • r ' 2 d*' 

(B.13) 

Taking the derivatives of Eqs. (B.12) and (B.13) in x' , 


K (d u e diO a e d 2 a e 

Id*- d*'J 2 (def 

(B.14) 

and 


K f du e du,^ a r d a a r 

Id* 1 d *'J" 2 (dx') a 

(B.15) 

Subtracting Eq. (B.9) from Eq. (B.7), 
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Multiplying Eq. (B.16) by K and substituting Eq. (B.14), 


E, 


(B.16) 


K 




=%^ +K: K-«.) Ar 


E, K) 2(d*') 


(B.17) 


Rearranging, 


d’o^ = 2K 


(dxT a e 


2s--2i-+(a e -a r )AT 


Solving Eq. (B.l) for a,, 


r *r 


(B.18) 


(B.19) 


Substituting Eq. (B.19) into Eq. (B.18), 

dV 


(dx 


Factoring and rearranging, 


^ = V’ " ff' + ft" + ( * AT 

) a e iK a r E r °r A r 


(B.20) 


d 2< *e OK 

fo r E r + «Al 

2X ' 

(dx-) a 

^ a r a c EJ£ e J 

°' = ®. 1 


Let the shear lag parameter £, which controls the thickness of the shear 
transfer region, be expressed as 


c 


I Ka e {cirE, + a c E e 
" 1 2 a r E r E c 


) 


Also let 




(B.22) 


(B.23) 


From the rule of mixtures, 


^E^drEr+O'E' 


(B.24) 
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and 


a 0 E 0 a 0 = a r E r a r + a c E c a c 


(B.25) 


Substituting Eqs. (B.24) and (B.25) into Eq. (B.23), X can be conveniently 
expressed as 


where <7^ , the far-field stress in the cracking layer, is 


^ = 


is. 

L$> 


<r a + E c (a # -a e )AT 


Substituting Eqs. (B.22) and (B.24) into Eq. (B.21) gives 
d*cr. 


(dx 


Placing Eq. (B.23) into Eq. (B.28), 


■f a! a. L a r E r ’> J 


ifSi-iSlo -X 

(d*')’ a S * 


This has a solution of the form, 


= Asinhf + Bcosh 

l <*« J 


' 2Qc' ) W c 

l * J V 

Applying boundary conditions <s e * 0 at x’- ±h , 

O.A8ud>(^) + Bcoah(^) + ^ 


0 = -Asinh 


K o. . 


+ £cosh 


(?) 


Xa\ 

4? 


Adding Eqs. (B.31) and (B.32) to find B, 


(B.26) 


(B.27) 


(B.28) 


(B.29) 


(B.30) 


(B.31) 


(B.32) 
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B = 


-Xo? 


4£*cosh(2£/i/a e ) 
Subtracting Eq. (B.32) from (B.31) and solving for A, 

A = 0 

Substituting Eqs. (B.33) and (B.34) into Eq. (B.30), 


o=- 


Xa 3 

C 


4f I coshf^' 

l<*. J 


cosh 


< a c > 


Xa 3 

4C 1 


Factoring, 


_Ao? 

4? 


cosh(2^c7a e ) 

cosh(2£7i/a e ) 


Substituting Eq. (B.26) into Eq. (B.36), 


<*c = 


cosh(2 £c’/a e ) 
cosh(2 CV a c) J 


Substituting Eq. (B.37) into Eq. (B.19), 

<r.a» o, 


<^r = 


cosh(2Cr7a c ) 
co8h(2£/i/a c ) J 


Or a , 

Substituting Eq. (B.37) into Eq. (B.7) and solving for 


o*. 

d u. 


du e _ 1 
d x'~E„ 


co8h(2C*yq e ) 

cosh(2£A/a c ) 


dx' 

a. + a AT 


(B.33) 


(B.34) 


(B.35) 


(B.36) 


(B.37) 


(B.38) 


(B.39) 


Placing Eq. (B.27) into Eq. (B.39), substituting Eqs. (B.24) and (B.25), 
factoring the loading terms, and integrating, 


a,E, [ a a g a J f'L coah( 2 frya e ) l , 

c " a 0 E 0 [a r E r + MT JJ 0 [ cosh( 2 £A/a e ) f* 

i r r i a ^ r ( 1 a<! 

J 0 [ a 0 E t { 2 £ cosh( 2 ^A/o e ) Jj 1 * 


+aAT 


(B.40) 


which gives 
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U . 3&J &&+,m.Ti a ‘ siph ( 2 ^'/ a ») 

• a„E„L“X ' I 2S.’cosh(2fA/a e )J 

, a A7V f 1 «*fl «■ 

* [ aX ( 2{r’ cosh(2fA/a.) J_ 

Substituting Eq. (B.38) into Eq. (B.9) and solving for 

OX 

+aA T (B.42) 

dx a r E r a^i cosh(2£A/a e ) J 

Placing Eq. (B.27) into Eq. (B.42), substituting Eqs. (B.24) and (B. 25), 
factoring the loading terms, and integrating, 


r a r E, [a n a a /„ _ wmll*-. cosh(2£r7®e)1 , _ A „,] J-t 

“'“I (it itt' " (0< ' '“^T" cosh(2fA/o«) J + “' AT r (B-43) 


(B.42) 


which gives 


o 0 <T a x' a c A? c fogO; / wyl ,_a t sinh(2£g /a e } + 

ar E r a 0 E 0 [a r E r K *'>“ }* 2£ cosh(2£A/a e ) J + <V “* 


' a,£ r a 0 E 0 [a, 
Factoring Eq. (B.44), 


(B.44) 


„ a„*' F F f. o. sinh(2{r'/g,) l 

r_ o,Ex[ ” £ * 2ft' cosh(2{A/o.) Jj 

&Tx T „ rv/ \( * a e 8inh(2{*7 a,.)^ 
+ ^[ a#£#ar + ae£ ^ ae " ar V"2Ct* cosh(2£A/a e ) J 


(B.45) 


The constants of integration in Eqs. (B.41) and (B.45) are set to zero to satisfy 
the displacement boundary conditions u r (0) = u e (0) = 0 . 


B.2 ENERGY EXPRESSIONS 

The energy criterion can be expressed as functions of the strain energy 
and external work alone in a Griffith energy balance: 
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(B.46) 


G = 


AW- ALT 


The condition for instantaneous appearance of a new crack is 

G £ G Ie (B.47) 

The change in internal strain energy, AXJ, has contributions from both the 
normal and shear stresses. The strain energy from normal stresses, U„ t is 


“•-fj 


od£ n dV 


(B.48) 


where 


E 

Substituting Eq. (B.49) into Eq. (B.48) and integrating over da, 


(B.49) 




E 


(B.50) 


The change in U„ when a new crack appears is 


(B.51) 


2U a \ k is the strain energy due to normal stresses of the volume element with 


cracks separated by a distance h. is the strain energy with cracks 

separated by 2h and can be found from 




|o, f ^*'+50. f 

2 J. E t 2 Jo 


(B.52) 


I II 

Note that due to symmetry, we can integrate from 0 to A and multiply the 
result by two, rather than integrating from -h to h. Substituting Eq. (B.38) 
into I, 
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b[if- 

f{ 


^£L + _?L_ 

2 a r E r 2a r E r 


1- 


2cosh 

a * , 

cosh 2 

r 2{*'Y| 

k J 

cosh 

r 2 Ch\ 

l J 

~ T 

cosh 2 

'2^1 
< a c J. 




g e Q a 

o r ^ r 


cosh 


2£x' 


1- 


cosh 


2£A 


cr.ff. 


(B.53) 


dx’ 


Substituting the trigonometric relation cosh 2 x = -^(cosh2x + 1) into Eq. (B.53) 

a 

and integrating, 


-a f* ^dx' = 
2 f Jo E r 




2a r E r 2a r E r 


A-^-tanhj^ 

C «e . 


'3 


sinh 


4&VA 
k«L > 2 


cosh 3 


q e q 0 

a r E r 


A-— ^tanh 

2C 


(¥1 


2fr) 

, a c J 




(B.54) 




Substituting Eq. (B.37) into II, 




2 E. 


2 cosh 


1- 


2£x 

CL. 


cosh 


ill coah-f^l 

i_J+ a '4 0 

? 


2# 

q 


Mx* (B.55) 


Substituting the trigonometric relation cosh 3 x = ^(cosh2x + 1) into Eq. (B.55) 
and integrating, 
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u r |dx= 
2 J. % 


2E„ 


A_i.tanh| 




sinh 


( 4£A > ' 

y a > 




0.J 


cosh 9 


2£A 

V fl c J 




Placing Eqs. (B.54) and (B.56) into Eq. (B.52), 

U I = £^ 
a? 




a r E c E r 


/i-^tanh 


2CA 

V fl c j 


-^-sinh 

'4£AV A 

8 C l 

. J 2 

cosh 2 

f 2 CM 
<*c J 


2a e a 0 <r a a e _ 


a r E r 


h S- tanh 

2C 


(sal 

V °c Jj 


(B.56) 


(B.57) 


The strain energy due to normal stresses in the volume element after a 
new crack forms and the cracks separated by a distance h is 


[ , f V* *2 - f V2 .2 

H t* ,+ H (B - 58) 

Due to symmetry, we integrate from 0 to hi 2 and multiply the result by two. 
Also, in Eq. (B.58), cr* and ff* are the normal stresses in the cracking layer 
and the rest of the laminate, respectively, derived with the boundary 
condition* a\ = 0 at x’- ±%. Then, by analogy with £7^, 
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u I = £ a l k 

" U 2aJ! r 


i °. g A 


a,.E,E r 


2a.o. 


'>.2, tanh f2* 

2 £ 1 - 


ijaU 

\ «£ L 0.1 « 




cosh 2 ! — 


a^r 


u 

— - tamh 


2 2£ 



*e J 


< 


0.0. 


Substituting Eqs. (B.57) and (B.59) into Eq. (B.51), 


(B.59) 


AC/ , = _!&&. 


a r®e^rC 


tanh 


(Ml 

-2tanh 


l O, J 


<°C J 


— sinh 
8 


4# 


N £A 1 • u 
+-2 — — sinh 

2a„ 4 




ii 

2a. 


coshl 2 ^' 


cosh 2 (^ 

l a c J 


UJ 






(B.60) 


Using the trigonometric identity sinh2x = 2sinhxcoshx, 



Simplifying, 


(B.61) 
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AU.-&&L 

' 2a r E r C 

— tanh] 

2 


2Ch 


— 3 tanh 


+4 


feMSM#-' 


( 2 ft} 

K a c ) 


-sech 2 


r C 


v°c y 




< (B.62) 


Placing Eq. (B.27) into Eq. (B.62), substituting Eqs. (B.24) and (B.25), and 
expanding the loading terms, 


V\ 


- [ <unh (f )- 2urfi (v)]} - 

Factoring, 

“••iaMHSl-MS'SWS)'"*'®) 

The strain energy contribution from shear stresses, U q> is given by 


(B.63) 


2 Ch 


\\ 




(B.64) 


o?E, 


“••111 


qdydV 


(B.65) 


where 


dy* 


dq 

G* 


(B.66) 


Substituting Eq. (B.66) into Eq. (B.65) and integrating over dq, 
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(B.67) 


v < 

y 

Gq. (B.67) is integrated over the volume of the shear transfer region, which 
has a total thickness of 2 a, since shear is transferred at both the top and 
bottom of the cracking layer. The change in U q when a new crack forms is 

At/,«[2L/,L-U,L] (B.68) 

2U q \ h is the strain energy due to shear stresses of the volume element with 
cracks separated by a distance h. U q \ , the strain energy with cracks 
separated by 2 h, can be found by solving Gq. (B.ll) for G**, and substituting 
the result into Gq. (B.67), giving 




2a q g 2 

Ka„ 


<&' 


(B.69) 


Simplifying, 




(B.70) 


Due to symmetry, we can once again integrate from 0 to h and multiply the 
result by two, rather than integrating from -h to h. Placing Gqs. (B.41) and 


(B.45) into Eq. (B.10) and squaring, 





sinh 2 J 

f 2&' 

“ 

\ J 



cosh 3 

f 2 CA" 





. 


o? 


(B.71) 


Substituting Gq. (B.71) into Gq. (B.70) and integrating, 
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uA = 


2Kal& 


,la * a 2 r E 2 c El 
Evaluating Eq. (B.75) from 0 to A, 


CM - &K 


a ( 

a] 4C 8inh > 

J 

-x’ 

cosh 1 

f 2£h 

°c > 

) 


ot 


(B.75) 


»» 4a iE^EfC 


~-sinh 


4C 


4 Ch 

a. 


?)-■ 


cosh 1 


:f2CM 

l J 


< 


(B.76) 


Using the trigonometric relation sinh2x = 2sinhxcoshx, Eq. (B.76) can be 
reduced to: 

KalajEl 


U. I. = 


[ 2 tanhf ^^l-^^8ech a f — llof (B.77) 

16 a?EfE?£* |_ l o t J o. I a, )\ ~ 

The strain energy due to normal stresses in the volume element when 
a new crack forms and with the cracks now separated by a distance h is 




(B.78) 


Due to symmetry, we integrate from 0 to h/2 and multiply the result by two. 
Also, in Eq. (4.33) q*is the shear stress derived with the boundary conditions 
a\ = 0 at x’ = ±%. Then, by analogy with U a ^ h , 


, 2XbX 
ajE*E* 


-5^-sinh 
a? 4f 1 

(*&} 
°c J 

-x’ 

8 ** cosh 2 

f g£ 

< a c J 





(B.79) 


Evaluating Eq. (B.79) from 0 to 
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| fo&K 4{ J 2 ^ 

,L_ 4 W coeh’f^l “ 


(B.80) 


Using the trigonometric relation sinh2x = 2sinhxcoshJC, Eq. (B.80) can be 


reduced to: 


17,1 = ^°;5 7 » l 2tanh f — V — aech»f oj. (B.81) 

* 1 6a.E,E r £ a e J a e l a «J 


,u \Qa 2 r E*EX L UJ <*« 

Substituting Eqs. (B.77) and (B.81) into Eq. (B.68), 


&U q = 


^ g-Lfnh fa' 
leo^EXC 3 [ U> 


-2tanhf^l^fsech»f^Vsech<^ 


(B.82) 

Solving Eq. (B.22) for If and placing into Eq. (B.82) and substituting Eqs. 
(B.27), (B.24), and (B.25), 


2tanh(^j-tanhg) + ^(s e ch*(^)-sech»(£)) 


(B.83) 


The total change in strain energy when a crack appears, AU, is simply 


AU = AU a + AU q 


(B.84) 


Substituting Eqs. (B.64) and (B.83) into Eq. (B.84), 

bXJ - ~ °C - a r) tqnhf j 


- ayM-^E]E r {a c -a f ) % Ar f J Wj)] (B . 8 5) 

2o r a.W L UJ l °c JJ 

The change in external work when a crack forms, AW, is given by 
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Aflr.[2n-wU ®»«> 

where W ^ and 2W[ A are the work done by the applied loading before and 
after the new crack forms, respectively. The former is found from 


= 2a 0 o a u r (h) 


(B.87) 


Evaluating Eq. (B.44) at x'=h, 


Ur{h)m ^.2& 

dfE, a 0 E ( 


^&L-(a e -a r )AT 


L a r&r 

Substituting Eq. (B.88) into Eq. (B.87), 

alalE'd - ala r a 0 E c E r ( a c - <0^(7. 

Za r a.E J E t 


h - ~ tanh 

to y 

JL 2 f 

< °e >. 


+ a r ATi (Bg8) 


2fr 


- tanh 


to 

l J] 


+ 2a 0 a r AToJx 


(B.89) 


The work done by the applied loading after the new crack forms is given by 

WV2a„<x„«;(|) (B.90) 

where is found from Eq. (4.22), derived using the boundary conditions 

<r* = 0 at x'= ±%. By analogy with W^, 


a c 2 a?g e qf - ala^E^ ( a e - a r )ATa a 
^a r a Q E r E 0 


a- t®nfa 

«c U;j 


(B.91) 


+ a 0 a r &Ta a h 


Substituting Eqs. (B.89) and (B.91) into Eq. (B.86), 


AW = 


2a]a]E e o 1 a - 2a i c a r E c E r a 0 {a e -a r ) q. A T 


2 Ca r a 0 E r E 0 


2 tanhf — j - tanh 

toY 

UeJ 

l °c J. 


(B.92) 
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The strain energy release rate is found by placing Eqs. (B.85) and (B.92) into 
Eq. (B.46), yielding, 


2tanh 


r 2C h\ 



< a c ). 


Placing Eqs. (B.24), (B.25), and (B.27) into Eq. (B.93), the strain 
release rate may also be expressed as 


o- 


2 tanh 


f- 

1 °./ 


-tanh 


(fl 


2 farErE' | 

Substituting Eq. (B.93) into Eq. (B.47), the cracking criterion becomes 




2 tanh 


-tohfsall 


< a c > 

l <*« JJ 


*G Ie 


(B.93) 

energy 

(B.94) 

(B.95) 
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APPENDIX C 


COMPUTER CODE MANUAL AND 
DOCUMENTATION 


4****>M**************************M***********««********************** 

USING CRACKOMATIC U 

CODE FOR THE PREDICTION OF 
MATRIX CRACKING 

c 1995 Jason R. Maddocks and Hugh L. McManus 
Massachusetts Institute of Technology 
Rm 33-311, 77 Massachusetts Ave. 

Cambridge MA 02139 (617) 253-0672 

*4^*******;MnM‘******4************************************************* 

Version 2.1 1/95 
WRITTEN IN MPW FORTRAN 
c 1988,1989 Language Systems Corp. 

***m^**************************************************************** 


This program is a research tool in the development stage 
and is supplied "as is" for the purpose of 
scientific collaboration. 

******** 4 ********************** 4 ************************************** 

L INTRODUCTION 

CRACKOMATIC II Version 2.1 calculates matrix crack density and reduced 
laminate properties in every ply of any arbitrary laminate as functions of 
temperature, thermal cycles, or any thermomechanical load history. 

n. INPUT FILES 

CRACKOMATIC II requires four kinds of prepared input files: 1) a laminate 
file containing material and layup information, 2) a fatigue toughness file 
for thermal cyclic loading, 3) a load profile for thermomechanical loading, and 
4) temperature/cycle dependent material property files (optional). These 
should be ASCII text files in the same folder or directory as the code. 
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Separate all entries by spaces and terminate every line, including the last 
one, with a carriage return. 

Ila. LAMINATE FILE FORMAT 

NPLY 

NMATi ANGLEi THICKi FLAGi 
repeat above line NPLY times 

EXi EYi NUXYi GXYi 

ALPHA li ALPHAS BETAli BETA2i 

repeat above 2 lines for each material 

NPLY is the number of ply groups in the laminate 
for each ply: 

NMATi is the material number of that ply group (number materials 
consecutively from 1) 

ANGLEi is the ply angle in degrees 
THICKi is the ply thickness 

FLAGi is a ply printing option. It should be set to 1 to print out the output 
data/results for ply t. Set to 0 if the printing of results for this ply is not 
needed. For each material: 

EXi is the longitudinal ply modulus 
EYi is the transverse ply modulus 
NUXYi is the mqjor Poisson's ratio 
GXYi is the shear modulus 
ALPHA li is the longitudinal ply CTE 
ALPHA2i is the transverse ply CTE 

BETAli and BETA2i are currently dormant ply CMEs- use 0.0 
EXAMPLE: (A P75/934 [G/45/90/-45]s laminate, English units) 

7 

1 .005 0 0 

1 .005 45 1 

1 .005 90 1 
1.010-451 

1.005 900 

1.006 45 0 
1.005 0 0 

34.3E6 0.9E6 .29 0.7E6 
-0.6E-6 16.E-6 0 0 

m>. FATIGUE TOUGHNESS FILE FORMAT 

NPOENTS 

NiGCi 

repeat above line NPOINTS times 
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NPOINTS is the number of data points that defines the Gc vs N curve for 
each point: 

Ni is the number of cycles 

GCi is the measured transverse fracture toughness corresponding to that 
number of cycles 

EXAMPLE: (P75/ERL1962, extrapolated from very limited data) 


6 


1 

.59566 

10 

.55627 

100 

.51947 

1000 

.48512 

10000 

.45303 

100000 

.42307 


lie. THERMOMECHANICAL LOAD PROFILE 

The thermomechanical load profile may accommodate any number of user- 
defined load increments, with thermal and mechanical loads applied either 
simultaneously or separately. Each increment is treated as a ramp, and 
loads are assumed to be quasi-static. Cycling effects are neglected. 

The code always starts at the stress free temperature and ramps down to the 
first load in the load profile, which accounts for cooldown at the end of the 
cure cycle. Thus the first entry in the user's load profile should be the 
ambient conditions, i.e. room temperature and zero mechanical stress. The 
initial cure cycle ramp is not included in the output. 

File format for thermomechanical loading: 

NPROF 

TiSIGMAi 

repeat above line NPROF times 

NPROF is the number of increments in the load profile. 

Ti and SIGMAi are the temperature and stress, respectively, at the i* 

increment, Ti is an actual temperature, not a AT. 

EXAMPLE: (English units) 

4 

70 0 

70 8.0E4 
-100 8.0E4 
-300 1.2E5 
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The laminate starts at ambient conditions, 70°F and zero applied mechanical 
stress. It is then loaded to 80 ksi with the temperature held constant. Next, 
with the mechanical stress held constant at 80 ksi, the temperature is 
reduced to -100°F. Finally, the temperature and mechanical stress are 
ramped simultaneously from the previous increment to -300°F and 120 ksi, 
respectively. 


Hd. TEMPERATURE/CYCLE DEPENDENT MATERIAL PROPERTY FILES 

The temperature dependent material file is similar to the fatigue toughness 
file format. All material constants are listed. 

File format: 

NPOINTS 

Ti EXi EYi NUXYi GXYi ALPHAli ALPHA2i GCi 
repeat above line NPOINTS times 

Ti is the temperature 

Other variables are as previously described 
EXAMPLE (P75/934) 

3 


-250 

.22791 

33.76E6 

.9E6 

.31 

1.1E6 

-0.43E-6 

21.923E-6 

75 

.22791 

41.97E6 

.83E6 

.35 

.61E6 

-0.584E-6 

19.18E-6 

250 

.22791 

45.36E6 

.81E6 

.30 

.46E6 

-0.365E-6 

26.455E-6 


File format for the material properties as function of cyclic loading is similar: 
NPOINTS 

Ni EXi EYi NUXYi GXYi ALPHAli ALPHA2i GCi 
repeat above line NPOINTS times 


m. INTERACTIVE SESSION 

Once you have defined some layups and material files, run CRACKOMATIC. 
The following capitalized text refers to the questions/options during the 
interactive session. 

The program will first ask for the laminate input file, give it one. Then it 
asks: 

COMPUTE MINIMUM (1), MAXIMUM (2), OR AVERAGE (3) CRACK DENSITY? 
Minimum seems to work best, so always choose it unless you specifically 
want to check theoretical maximum or average densities. 
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The code then outputs a review of the laminate, and a complete set of 
calculated laminate properties. 

ANALYSIS TYPE- 

1 = CRACK DENSITY AND PROPERTIES AS FUNCTION OF DELTA-T 

2 = CRACK DENSITY AND PROPERTIES AS FUNCTION OF N 

3 = CRACK DENSITY FROM THERMOMECHANICAL LOAD PROFILE 

4 = READ NEW LAMINATE 

5 = QUIT 

Choices 1, 2, and 3 have their own sections below. 

Choice 4 lets you pick a new laminate and/or change your choice of crack 
density. Choice 5 stops the code, leaving the session in an editable text 
window, where the results can be cut and pasted into other documents. 


Ilia. ANALYSIS TYPE 1 

This analysis calculates the progressive change in crack density and laminate 
properties as functions of decreasing temperature. The code asks: 

GIVE G (TRANSVERSE PLY FRACTURE TOUGHNESS), 

SHEAR LAG FACTOR 

AND LAMINATE STRESS FREE TEMPERATURE 
G is the critical strain energy release rate. 

SHEAR LAG FACTOR is a geometric parameter which can reasonably range 
from around 0.5 to around 2.0. 

LAMINATE STRESS FREE TEMPERATURE is usually the laminate cure 
temperature. The code then asks 

GIVE TEMPERATURE RANGE AND INCREMENT: 

INITIAL TEMP, FINAL TEMP, AND TEMP INCREMENT 
User option variables to control the printed output 

WANT TO SOFTEN LAMINATE AS IT PROGRESSES? 

User option to include material softening effect (answer yes or no) 

WANT TO INCORPORATE TEMPERATURE DEPENDENCE OF MATERIAL 
CONSTANTS (REQUIRES INPUT FILE)? 

User option to include temperature dependent material properties, required 
input file with appropriate data (see format in Section lid.) 

The output is a tab-separated table of temperatures, crack densities and 
laminate lon gi t udinal stiffness and CTE. These can be used to generate plots 
of progressive cracking and changin g l amin ate properties as the temperature 
is progressively decreased. At the final temperature, the program computes 
all the degraded laminate properties. 
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nib. ANALYSIS TYPE 2 


This analysis calculates the change in crack density and laminate properties 
as functions of numbers of constant thermal cycles. The code asks: 

INPUT GKN) FILE NAME, OR HELP IP YOU NEED IT 

Give it the name of a FATIGUE TOUGHNESS FILE (see format in Section 
Hb.). 

GIVE SHEAR LAG FACTOR AND GREATEST DELTA-T 

Give the shear lag factor described above, and the DELTA-T at the lowest 
temperature in the cycle. (The highest temperature or cycle R-value are 
assumed to be the same as those used to generate the data in the fatigue 
toughness file, so you don't input either of these). 

INCREASE N LINEARLY (ENTER Y) OR EXPONENTIALLY (N)? 

This choice controls the output. A "Y” will give output suitable for making a 
linear plot, while an "N" will generate output suitable for making a semi-log 
plot. If you choose Y, you are asked 

GIVE MAXIMUM N AND INCREMENT 

which is self-explanatory; if you choose N, you are asked 

GIVE MAXIMUM N AND POINTS PER DECADE 

which is almost so; points per decade is the number of plot 

points generated for each power often on the plot. 

WANT TO SOFTEN LAMINATE? 

User option to include material softening effect 

WANT TO INCORPORATE MATERIAL CONSTANTS AS FUNCTION OF THERMAL 
CYCLES (REQUIRES INPUT FILE)? 

User option to include material properties as a function of thermal cycles, 
required input file with appropriate data. 

The output is a tab-separated table of number of cycles, crack densities, and 
laminate longitudinal stiffness and CTE. These can be used to generate plots 
of cracking and chang ing l amina te properties as the laminate is thermally 
cycled. 

me. ANALYSIS TYPE 3 

This analysis calculates the progressive change in crack density and laminate 
properties due to a user-defined thermomechanical load history. The code 
first asks: 

PLEASE INPUT LOAD HISTORY FILE 

Enter the name of the THERMOMECHANICAL LOAD PROFILE file (see 
format in Section He.) 

GIVE G (TRANSVERSE PLY FRACTURE TOUGHNESS), 
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SHEAR LAG FACTOR 

AND LAMINATE STRESS FREE TEMPERATURE 
G is the critical strain energy release rate. 

SHEAR LAG FACTOR is a geometric parameter which can reasonably range 
from around 0.5 to around 2.0. 

LAMINATE STRESS FREE TEMPERATURE is usually the laminate cure 
temperature. The code then asks 

WANT TO SOFTEN LAMINATE AS IT PROGRESSES? 

User option to include material softening effect 

WANT TO INCORPORATE TEMPERATURE DEPENDENCE OF MATERIAL 
CONSTANTS (REQUIRES INPUT FILE)? 

User option to include temperature dependent material properties, required 
input file with appropriate data (see format in Section Ed.) 

The output is a tab-separated table of temperatures, applied mechanical 
stresses, crack densities and laminate longitudinal stiffness and CTE. These 
can be used to generate plots of progressive cracking and changing laminate 
properties with progressive changes in temperature, stress, or both. At the 
final load, the program computes all the degraded laminate properties. 


IV. A FEW USEFUL TRICKS 

In the Macintosh version of the code, the output remains in an editable, 
saveable text window following execution. The edit and file commands are 
functional. Usually, all work in a session can be accessed by scrolling; very 
long sessions save themselves in a file and must be closed and reopened using 
an editor. 

The outputs of Analysis Types 1, 2, and 3 are in tab-separated tables that can 
be cut and pasted directly into spreadsheets or plotting programs. The code 
expects all inputs in a consistent set of units. Make sure, for example, that if 
you are using metric units you enter ply thicknesses in METERS, not 
millimeters as they are usually reported. 

A "hack" in the code allows you to give all properties in English units except 
G fc (which is usually reported in metric) by entering negative G tc values; the 
code converts the absolute value of G k to English units. 


V. SAMPLE SESSION 

Three samples sessions are listed below. The first is an example of Analysis 
Type 1. This example includes material softening effects and temperature 
dependent material properties. The second session is an example of Analysis 
Type 2 using material softening effects. The third shows an Analysis Type 3 
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including material softening effects and temperature dependent material 
properties. User input is in bold italics.. 

Va. SAMPLE SESSION A: ANALYSIS TYPE lQE_ET5/934 r0/45/90/-451a 


INPUT LAYUP PILE OR HELP IF YOU NEED IT 
p75934*fltf 


PLY 

LAMINATE 

MATERIAL THICKNESS 

ANGLE 

1 

1 

0.0060 

0.00 

2 

1 

0.0060 

45.00 

3 

1 

0.0050 

90.00 

4 

1 

0.0100 

•45.00 

5 

1 

0.0060 

90.00 

6 

1 

0.0060 

46.00 

7 

1 

0.0060 

0.00 


MATERIAL PROPERTIES FOR MATERIAL 1- 
Ell 3.4300E+07 E22 9.0000E+05 

NU12 2.9000E-01 G12 7.0000E+05 

ALPHA1 -6.0000E-07 ALPHA2 1.6000E-06 
BETA1 O.OOOOE+OO BETA2 0.0000E+00 

ENGINEERING CONSTANTS FOR UNCRACKED LAMINATE- 
Ella 1.2317E+07 E22a 1.2317E+07 G a 4.6943E+06 
V12 a 3.U91E-01 V21a 3.1191E-01 

ALPHA1 a -6.0484E-08 ALPHA2 a -6.0484E-08 ALPHA12 a -4.8557E-19 
BETA1 a O.OOOOE+OO BETA 2 a O.OOOOE+OO BETA12 = 0.0000E+00 

COMPUTE: (1) MINIMUM, (2) MAXIMUM, (3) AVERAGE CRACK SPACING 


ANALYSIS TYPE- 

1 a CRACK DENSITY AS FUNCTION OF DELTA-T 

2 a CRACK DENSITY AS FUNCTION OF CYCLES (N) 

3 « CRACK DENSITY FROM THERMOMECHANICAL LOAD PROFILE 

4 a READ NEW LAMINATE 

5 a QUIT 
1 


GIVE G (TRANSVERSE PLY FRACTURE TOUGHNESS) 
SHEAR LAG FACTOR, 

LAMINATE STRESS FREE TEMPERATURE 
-40.65350 

GIVE TEMPERATURE RANGE AND INCREMENT: 
INITIAL TEMP, FINAL TEMP, AND TEMP INCREMENT 
75-25025 

WANT TO SOFTEN LAMINATE AS IT PROGRESSES? 

y 
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WANT TO INCORPORATE TEMPERATURE DEPENDENCE OF 
MATERIAL CONSTANTS (REQUIRES INPUT FILE) ? 
y 


PLEASE INPUT TEMPERATURE DEPENDENT MATERIAL FILE 
p75934_tdep 

TRANSVERSE PLY FRACTURE TOUGHNESS = 4.0000E+01 (METRIC) 
TRANSVERSE PLY FRACTURE TOUGHNESS = 2.279 IE-01 
SHEAR LAG FACTOR* 0.65 
STRESS FREE TEMPERATURE * 350.00 



* PLY 2 * 

* PLY 3 * 

•PLY 4* 

LAMINATE 

Temperature 

Crack Density 

Crack Density 

Crack Density 

Stiffness 

CTE 

75 

0.0000E+00 

O.OOOOE+OO 

O.OOOOE+OO 

1.4776E+07 

-7.3510E-08 

50 

0.0000E+00 

O.OOOOE+OO 

8.0141E+00 

1.4571E+07 

-5.7974E-08 

25 

O.OOOOE+OO 

O.OOOOE+OO 

1.5183E+01 

1.4367E+07 

-4.1396E-08 

0 

O.OOOOE+OO 

O.OOOOE+OO 

1.8989E+01 

1.4173E+07 

-1.9946E-08 

-25 

O.OOOOE+OO 

O.OOOOE+OO 

2.2130E+01 

1.398 1E+07 

2.6973E-09 

-50 

O.OOOOE+OO 

O.OOOOE+OO 

2.4945E+01 

1.3788E+07 

2.6106E-08 

-75 

2.6459E+01 

3.6858E+01 

2.7567E+01 

1.3517E+07 

-1.339 IE-08 

-100 

3.3482E+01 

4.6949E+01 

3.0048E+01 

1.3300E+07 

-8.4506E-09 

-125 

3.9021E+01 

5.4842E+01 

3.2430E+01 

1.3086E+07 

-1.4497E-10 

-150 

4.3892E+01 

6.1767E+01 

3.4737E+01 

1.2872E+07 

9.5830E-09 

-175 

4.837 1E+01 

6.8126E+01 

3.6984E+01 

1.2660E+07 

2.0111E-08 

-200 

5.2588E+01 

7.4108E+01 

3.9184E+01 

1.2446E+07 

3.1153E-08 

-225 

5.6616E+01 

7.982 1E+01 

4.1345E+01 

1.2233E+07 

4.2564E-08 

-250 

6.0469E+01 

8.5284E+01 

4.3457E+01 

1.2027E+07 

5.3318E-08 


ALL LAMINATE PROPERTIES AT FINAL TEMPERATURE 
Ell* 1.2027E+07 E22* 1.2026E+07 G * 4.628 1E+06 
V12 * 3.0888E-01 V21 * 3.0887E-01 

ALPHA1 * 5.3318E-08 ALPHA2 * 5.255 IE-08 ALPHA12 = -4.8800E-08 
BETA1 * 0.0000E+00 BETA2 * O.OOOOE+OO BETA12 * 0.0000E+00 

ANALYSIS TYPE- 

1 = CRACK DENSITY AS FUNCTION OF DELTA-T 

2 » CRACK DENSITY AS FUNCTION OF CYCLES (N) 

3 * CRACK DENSITY FROM THERMOMECHANICAL LOAD PROFILE 

4 * READ NEW LAMINATE 
5* QUIT 

5 


Vb SAMPLE SESSION B : ANALYSIS TYPE 2..QF P75/EBL196 2 

r0/90/0/901. 

INPUT LAYUP FILE OR HELP IF YOU NEED IT 
p75*nzn 
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PLY 

MATERIAL 

LAMINATE 

THICKNESS 

ANGLE 

1 

1 

0.0050 

0.00 

2 

1 

0.0050 

90.00 

3 

1 

0.0050 

0.00 

4 

1 

0.0100 

90.00 

5 

1 

0.0050 

0.00 

6 

1 

0.0050 

90.00 

7 

1 

0.0050 

0.00 


MATERIAL PROPERTIES FOR MATERIAL 1- 
Ell 3.4300E+07 E22 9.0000E+05 

NU12 2.9000E-01 G12 7.0000E+05 

ALPHA1 -5.3000E-07 ALPHA2 2.2000E-05 
BETA1 O.OOOOE+OO BETA2 O.OOOOE+OO 

ENGINEERING CONSTANTS FOR UNCRACKED LAMINATE- 
Ell- 1.7635E+07 E22 a 1.7635E+07 G * 7.0000E+05 
V12* 1.4830E-02 V21 = 1.4830E-02 

ALPHA1 * 2.0225E-07 ALPHA2 a 2.0225E-07 ALPHA12 a -6.8180E-18 
BETA1 a O.OOOOE+OO BETA2 = O.OOOOE+OO BETA12 a O.OOOOE+OO 

COMPUTE: (1) MINIMUM, (2) MAXIMUM, (3) AVERAGE CRACK SPACING 
1 

ANALYSIS TYPE- 

1 a CRACK DENSITY AS FUNCTION OF DELTA-T 

2 a CRACK DENSITY AS FUNCTION OF CYCLES (N) 

3 a CRACK DENSITY FROM THERMOMECHANICAL LOAD PROFILE 

4 a READ NEW LAMINATE 

5 a QUIT 

2 


INPUT G(N) FILE NAME OR HELP IF YOU NEED IT 
p75g* 

GIVE SHEAR LAG FACTOR AND GREATEST DELTA-T 
.65-600 

INCREASE N LINEARLY (ENTER Y) OR EXPONENTIALLY (N)? 

y 


GIVE MAXIMUM N AND INCREMENT 
50025 

WANT TO SOFTEN LAMINATE? 

y 


WANT TO INCORPORATE MATERIAL CONSTANTS 
AS FUNCTION OF CYCLES? 
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It 


FRACTURE TOUGHNESS AS A FUNCTION OF CYCLE 
N FRACTURE TOUGHNESS 
1.0000E+00 5.9566 E-01 

1.0000E+01 5.5627E-01 

1.0000E+02 5. 1947E-01 

1.0000E+03 4.8512E-01 

1.0000E+04 4.5303 E-01 

1.0000E+05 4.2307E-01 

SHEAR LAG FACTOR^ 0.65 
MAXIMUM DELTA-T s -600.00 



* PLY 1 * 

* PLY 2 * 

* PLY 4 * 

LAMINATE 

Cycle 

Crack Density 

Crack Density 

Crack Density 

Stiffness 

CTE 

0 

1.8201E+01 

1.8037E+01 

3.5615E+01 

1.7558E+07 

9.588 IE-08 

25 

3.4207E+01 

3.4483E+01 

3.7763E+01 

1.7540E+07 

6.3889E-08 

50 

3.5999E+01 

3.6245E+01 

3.8244E+01 

1.7538E+07 

5.9886E-08 

75 

3.7677E+01 

3.7901E+01 

3.8733E+01 

1.7536E+07 

5.6050E-08 

100 

3.9273E+01 

3.9480E+01 

3.923 1E+01 

1.7533E+07 

5.2329E-08 

125 

3.9423E+01 

3.9629E+01 

3.9280E+01 

1.7533E+07 

5.1975E-08 

150 

3.9568E+01 

3.9773E+01 

3.9326E+01 

1.7533E+07 

5.1633E-08 

175 

3.9712E+01 

3.9916E+01 

3.9374E+01 

1.7533E+07 

5.1291E-08 

200 

3.9856E+01 

4.0059E+01 

3.942 1E+01 

1.7533E+07 

5.0950E-08 

225 

4.0000E+01 

4.0201E+01 

3.9468E+01 

1.7532E+07 

5.0609E-08 

250 

4.0143 E+01 

4.0343E+01 

3.9516E+01 

1.7532E+07 

5.0269E-08 

275 

4.0286E+01 

4.0485E+01 

3.9563E+01 

1.7532E+07 

4.9930E-08 

300 

4.0429E+01 

4.0626E+01 

3.9610E+01 

1.7532E+07 

4.9591E-08 

325 

4.0572E+01 

4.0767E+01 

3.9658E+01 

1.7532E+07 

4.9252E-08 

350 

4.0713E+01 

4.0907E+01 

3.9706E+01 

1.753 1E+07 

4.8915E-08 

375 

4.0854E+01 

4.1047E+01 

3.9754E+01 

1.753 1E+07 

4.8575E-08 

400 

4.0995E+01 

4.1187E+01 

3.9802E+01 

1.7531E+07 

4.8239E-08 

425 

4.1135E+01 

4. 1327 E+01 

3.9850E+01 

1.753 1E+07 

4.7902E-08 

450 

4.1275E+01 

4.1465E+01 

3.9898E+01 

1.753 1E+07 

4.7566E-08 

475 

4.1415E+01 

4.1604E+01 

3.9946E+01 

1.7530E+07 

4.7230E-08 

500 

4.1555E+01 

4.1742E+01 

3.9995E+01 

1.7530E+07 

4.6894E-08 


ALL LAMINATE PROPERTIES AT FINAL TEMPERATURE 
Ells L7530E+07 E22 = 1.7564E+07 G = 5.6136E+05 
V12 s 1.1942E-02 V21= 1.1965E-02 

ALPHA1 s 4.6894E-08 ALPHA2 * 7.6861E-08 ALPHA12= -7.6337E-18 
BETA1 s O.OOOOE+OO BETA2 = O.OOOOE+OO BETA12 * O.OOOOE+OO 

ANALYSIS TYPE- 

1 « CRACK DENSITY AS FUNCTION OF DELTA-T 

2 * CRACK DENSITY AS FUNCTION OF CYCLES (N) 

3 s CRACK DENSITY FROM THERMOMECHANICAL LOAD PROFILE 

4 s READ NEW LAMINATE 

5 = QUIT 
5 
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input layup file or help if you need it 

p75934xftif 


PLY 

LAMINATE 

MATERIAL THICKNESS 

ANGLE 

1 

1 

0.0060 

0.00 

2 

1 

0.0060 

45.00 

3 

1 

0.0060 

90.00 

4 

1 

0.0100 

-45.00 

5 

1 

0.0060 

90.00 

6 

1 

0.0050 

46.00 

7 

1 

0.0050 

0.00 


MATERIAL PROPERTIES FOR MATERIAL 1- 
Ell 3.4300E+07 E22 9.0000E+05 

NU12 2.9000E-01 G12 7.0000E+05 

ALPHA1 -6.0000E-07 ALPHA2 1.6000E-06 
BETA1 0.0000E+00 BETA2 O.OOOOE+OO 

ENGINEERING CONSTANTS FOR UNCRACKED LAMINATE- 
Ell s 1.2317E+07 E22 = 1.2317E+07 G = 4.6943E+06 
V12 * 3.119 IE-01 V21 * 3.1191E-01 

ALPHA1 - -6.0484E-08 ALPHA2 a -6.0484E-08 ALPHA12 = -4.8657E-19 
BETA1 a O.OOOOE+OO BETA2 a 0.0000E+00 BETA 12 a O.OOOOE+OO 

COMPUTE: (1) MINIMUM, (2) MAXIMUM, (3) AVERAGE CRACK SPACING 


ANALYSIS TYPE- 

1 = CRACK DENSITY AS FUNCTION OF DELTA-T 

2 = CRACK DENSITY AS FUNCTION OF CYCLES (N) 

3 = CRACK DENSITY FROM THERMOMECHANICAL LOAD PROFILE 

4 a READ NEW LAMINATE 

5 a QUIT 
3 


PLEASE INPUT LOAD HISTORY FILE 
p75984_pro(U» 

GIVE G (TRANSVERSE PLY FRACTURE TOUGHNESS) 

SHEAR LAG FACTOR, 

LAMINATE STRESS FREE TEMPERATURE 
-40.65350 

WANT TO SOFTEN LAMINATE AS IT PROGRESSES? 

y 

WANT TO INCORPORATE TEMPERATURE DEPENDENCE OP 
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MATERIAL CONSTANTS (REQUIRES INPUT PILE) ? 

y 


PLEASE INPUT TEMPERATURE DEPENDENT MATERIAL FILE 
p75934Jdep 

TRANSVERSE PLY FRACTURE TOUGHNESS a 4.0000E+01 (METRIC) 
TRANSVERSE PLY FRACTURE TOUGHNESS = 2.2791E-01 
SHEAR LAG FACTOR * 0.66 

STRESS FREE TEMPERATURE = 350.00 


•PLY 2* * PLY 3 • * PLY 4 • LAMINATE 

Temp Applied Stress Crack Density Crack Density Crack Density Stiffness CTE 


75 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

1.4776E+07 

-7.3510E-08 

50 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

8.014E+00 

1.4671E+07 

-5.7974E-08 

25 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

1.518E+01 

1.4367E+07 

-4.1396E-08 

0 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

1.899E+01 

1.4173E+07 

-1.9946E-08 

-25 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

2.213E+01 

1.398 1E+07 

2.6973E-09 

-50 

O.OOOE+OO 

O.OOOE+OO 

O.OOOE+OO 

2.495E+01 

1.3788E+07 

2.6106E-08 

-75 

O.OOOE+OO 

2.646E+01 

3.686E+01 

2.757E+01 

1.3517E+07 

-1.3391E-08 

-100 

O.OOOE+OO 

3.348E+01 

4.695E+01 

3.005E+01 

1.3300E+07 

-8.4506E-09 

•125 

1.143E+04 

4.319E+01 

6.3UE+01 

3.438E+01 

1.3061E+07 

-1.6840E-08 

•150 

2.286E+04 

5.129E+01 

7.633E+01 

3.847E+01 

1.2828E+07 

-2.1044E-08 

-175 

3.429E+04 

5.862E+01 

8.825E+01 

4S40E+01 

1.2596E+07 

-2.3427E-08 

•200 

4.571E+04 

6.549E+01 

9.940E+01 

4.622E+01 

1.2365E+07 

-2.4731E-08 

•225 

5.714E+04 

7.207E+01 

1.101E+02 

4.996E+01 

1.2134E+07 

•2.5275E-08 

•250 

6.857E+04 

7.840E+01 

1S04E+02 

5.360E+01 

1.1912E+07 

-2.6036E-08 

•275 

8.000E+04 

8.357E+01 

1S88E+02 

5.661E+01 

1.1883E+07 

-4.5323E-08 


ALL LAMINATE PROPERTIES AT FINAL TEMPERATURE 
Ell a 1.1883E+07 E22* 1.1910E+07 G a 4.5773E+06 
V12= 3.1317E-01 V21 a 3.1389E-01 

ALPHAl a -4.5323E-08 ALPHAS a 4.9401E-09 ALPHA 12 = -4.7435E-08 
BETA1 a 0.0000E+00 BETA2 a O.OOOOE+OO BETA12 a O.OOOOE+OO 


ANALYSIS TYPE- 

1 a CRACK DENSITY AS FUNCTION OF DELTA-T 

2 a CRACK DENSITY AS FUNCTION OF CYCLES (N) 

3 a CRACK DENSITY FROM THERMOMECHANICAL LOAD PROFILE 

4 « READ NEW LAMINATE 
6 a QUIT 

5 
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APPENDIX D 


CRACK DENSITY DATA AND ANALYTICAL 

PREDICTIONS 


This appendix presents all results of the progressive thermal and 
mechanical loading tests of the [O 4 / 454 / 9 O /45J„ [0 a /45a/90a/-45J„ and 
[0 a /60 a /-60 2 ] i laminates. Thermal results comprise crack data from both 
interior and edge inspections. Average observed crack density is plotted as a 
function of applied load or decreasing temperature. Error bars represent one 
standard deviation. Predicted crack density curves are plotted on the same 
graphs as the experimental data. 
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[0/45/90/-45 1 

4 4 4 4 s 



Figure D.l Experimental results and analytical predictions of crack 
density vs. progressive applied load. 45 4 ply group of 
[O 4 / 454 / 9 O 4 /- 454 L laminate. 
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[0 745 7907-45 1 

4 4 4 4 s 

45, Ply Group 

4 



Temperature (°C) 


Figure D.2 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 45* ply 
group of [OV 454 / 9 O 4/-45 J, edge specimens. 


187 



[0/45/907-45 1 

4 4 4 4 8 



Temperature (°C) 


Figure D.3 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 45« ply 
group of[0V45V904/-45J. interior specimens. 
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[0/45/907-45 1 

4 4 4 4 s 



Figure D.4 Experimental results and analytical predictions of crack 
density vs. progressive applied load. 90< ply group of 
[O 4 / 45 / 9 O 4 M 5 J, laminate. 
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[07457907-45 1 

4 4 4 4 S 



Temperature (°C) 


Figure D.5 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 90* ply 
group of[0y454/90/45J. edge specimens. 
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[0/45/90/-45 1 

4 4 4 4 S 



Temperature (°C) 


Figure D.6 Experimental results and analytical predictions 
density vs. progressively decreasing temperature, 
group of[0V454/904/-45J. interior specimens. 


of crack 
90 4 ply 
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[07457907-45,] 
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Figure D.7 Experimental results and analytical predictions of crack 
density vs. progressive applied load. -45* ply group oi 
[O 4 / 454 / 9 O 4 /- 45 J, laminate. 
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[0/45/90 /-45 1 

4 4 4 4 s 



Temperature (°C) 


Figure D.8 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. -45* ply 
group of [ 0 ^ 45 ^ 90 ^- 454 ], edge specimens. 
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[0/45/90/-45J 

4 4 4 4 s 
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Figure D.9 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. -45$ ply 
group of [Oy454/904/-45 J, interior specimens. 
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Figure D.10 Experimental results and analytical predictions of crack 
density vs. progressive applied load. 45 a ply group of 
[Oa/45a/90a/-452], laminate. 
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Figure D.ll Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 45* ply 
group of [0a/45^90a/-45j. edge specimens. 
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Figure D.12 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 45* ply 
group of[0*/45*/9<y-45*]. interior specimens. 
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Figure D.13 Experimental results and analytical predictions of crack 
density vs. progressive applied load. 90* ply group of 
[0*/45*/90^-45*]. laminate. 
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Figure D.14 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 90 a ply 
group of [0^45^90 a/-45J, edge specimens. 
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Temperature (°C) 


Figure D.15 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 90* ply 
group of [<y45a/90a/-45d. interior specimens. 
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0 10 20 30 40 50 


Applied Load (kN) 


Figure D.16 Experimental results and analytical predictions of crack 
density vs. progressive applied load. -454 ply group of 
[0 3 /45^90 9 /-45 2 ]. laminate. 
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[0 2 /45 2 /90 2 /-45 2 ] 8 
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Figure D.17 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. -404 ply 
group of [03/453/90/45*]. edge specimens. 
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Figure D.18 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. -45 4 ply 
group of [0^45^90a/-45J. interior specimens. 
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Figure D.19 Experimental results and analytical predictions of crack 
density vs. progressive applied load. 60 a ply group of [(V60a/- 
60 J. laminate. 
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Figure D.20 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 60 a ply 
group of[0^W-60 a ], edge specimens. 


I 
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Figure D.21 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. 60* ply 
group of[0^W-60J. interior specimens. 
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[ 0 2 / 60 /- 60 2 ] s 



Figure D.22 


Experi me ntal results and analytical predictions of crack 
density vs. progressive applied load. - 60 * ply group of 
[ 0 a/ 60 a/- 60 J. laminate. 
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Figure D.23 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. *604 ply 
group of [0^60y-60J. edge specimens. 
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Figure D .24 Experimental results and analytical predictions of crack 
density vs. progressively decreasing temperature. -0O4 ply 
group of[<V 6 <y- 60 J. interior specimens. 
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APPENDIX E 


FIRST-PLY-FAILURE STRESS-STRAIN DATA 


This appendix presents stress-strain data for the [0/54/90/45*],, 
[09/459/90/453],, and [Og/60/603], tensile coupons. Stress-strain curves for 
the loadings after which first ply failure was observed are plotted. The 
specimen name is included with each plot. The first character indicates the 
layup: "E", "F", and "G" represent [O3/6O/6O3]., [03/453/90/453],, and 

[04/454/90/45*],, respectively. The second character represents the panel 
from which the specimen was cut, and the third indicates the specimen's 
position on the cutting plan in Figure 5.1. The final character is either ”N", 
"W", or "M M : narrow thermal specimen, wide thermal specimen, or tensile 
coupon, respectively. 
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Stress (MPa) 


Specimen GA-1M 



Figure E.l Stress-strain plot for GA-1M [O4/454/9O4M5J, tensile coupon 
progressively loaded to 15.6 kN. 
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Specimen GA-2M 



Figure E.2 Stress-strain plot for GA-2M [(y45*/9(V-45J. tensile coupon 
progressively loaded to 15.6 kN. 
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Stress (MPa) 


Specimen GA-3M 



Figure E.3 Stress-strain plot for GA-3M [0,/45 < /90 ) /-45J. tensile coupon 
progressively loaded to 15.6 kN. 
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Specimen GB-2M 



Strain (jxe) 


Figure E.4 Stress-strain plot for GB-2M tensile coupon 

progressively loaded to 15.6 kN. 
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Stress (MPa) 




Specimen FA-2M 



Strain (pe) 


Figure E.6 Stress-strain plot for FA-2M [0a/45a/90a/-45a], tensile coupon 
progressively loaded to 17.8 kN. 


216 



Stress (MPa) 


Specimen FA-3M 



Strain (pe) 


Figure E.7 Stress-strain plot for FA-3M i [0 a /45 a /90 a /-45aL tensile coupon 
progressively loaded to 17 .8 kN. 



Specimen FB-1M 



Strain (joe) 


Figure E.8 Stress-strain plot for FB-1M [0a/45j/90a/-45J, tensile coupon 
progressively loaded to 13.3 kN. 
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Specimen FB-2M 



Strain (jxe) 


Figure E.9 Stress-strain plot for FB-2M [0^45a/90a/-45 2 ]. tensile coupon 
progressively loaded to 15.6 kN. 
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Specimen FB-3M 
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Figure E.10 Stress-strain plot for FB-3M [0a/45a/90a/-45j], tensile 
progressively loaded to 17.8 kN. 


coupon 
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Figure E.12 Stress-strain plot for EA-2M [0a/60a/-60d« tensile coupon 
progressively loaded to 15.6 kN. 
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Stress (MPa) 





Specimen EB-1M 



Strain (ji£) 


Figure E.14 Stress-strain plot for EB-1M [(V6(V-60i], tensile coupon 
progressively loaded to 13.3 kN. 
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